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ABSTRACT 


Composite structure repair methodology has been developed to 
specific applications (typically in small area and limited to secondary 
structure) and is being extended to Large Area Composite Structure 
Repair (with target extension to primary structures). Therefore, the 
repair becomes more critical because we get redistribution of stresses 
that can also affect the zones outside of the repair area. For this reason, an 
analytic evaluation of the repair's reliability has to be performed to 
define a parameter which reflects on the effectiveness of the repair. 

In this work, we establish a principal guideline to evaluate the 
redundancy and compare the reliability of the repair to the reliability 
of the parent structure (i.e., the structure in the undamaged state). 
The approach adopted is to utilize structural finite element analysis to 
compute the state of stress at all the spatial elements of structure of 
the damaged state and of the candidate repaired state. The reliability 
of these two spatially non-uniform stresses is computed by a probabil- 
istic failure criterion. Thus, we can optimize the repair configuration by 
varving the strength and the stiffness of any element in the repair site 


by varying lamination angles, and selectively using hybrid materials. 
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I. INTRODUCTION 


A. BACKGROUND 

Composite materials are being utilized for commercial and mili- 
tary applications. Their combination of superior mechanical properties 
(especially the redundancy), design versatility, and low mass density 
make them attractive structural materials. Principal applications of 
composites have been on aircraft, especially on fixed-wing aircraft 
[So F-18, F-16, F-15, F-14). 

As the structural use of composite material on aircraft increases, 
procedures to repair damage become essential. Composite repair pro- 
grams have been undertaken both by commercial and military sources 
1 Repair concepts have been developed which include bolted 
patches repairs and bonded patches repairs. At the beginning. the 
repair technologv was targeted for specific applications, typically lim- 
ceto small areas in secondary structures. More recently, repair 
techniques have been extended to larger areas, with target extension 
to primary structures (which inherently have high strength and dura- 
bilitv requirements). The particular techniques developed to repair a 
given component are dependent on the availability of specific facilities 
and the level of structural integrity. 

ict Reference 1) repair procedures were developed for two 
repair configurations. The configurations that were considered are a 


scarf inclusion and an external patch. A portable repair kit was 


developed to provide vacuum pressure and heat in an integral rubber 
blanket. These techniques permit repairs, verified by testing, that are 
capable of restoring between 80 percent and 100 percent of the par- 


ent laminate ultimate allowable loading. 


B. STRUCTURAL REPAIR DESIGN 

Repair methodology development through testing has several 
inherent limitations. Only a specific structural configuration can be 
examined and optimized for variation of repair configuration in terms 
of the resultant effectiveness. A large number of experiments are fre- 
quentlv impractical due to economic and time constraints. This calls 
for the establishment of an analytical principle to guide repair 
optimization. 

In designing repair of composite structures, the external loading 
is known and the objective is to restore the structural functionality 
while recovering the overall structural reliability, which includes both 
the repair and the parent structure. The stiffness of the repair See 
ment and the strength of its attachment may both influence the overall 
Structural reliability by altering the state of stresses in the parent 
structure. Thus, a repair segment with high stiffness and strong 
attachment may cause structural failure elsewhere, resulting in a 
strong repair but a weaker overall structure. On the other hand, a 
repair with proper stiffness but with insufficient attachment strength 
may fail in the attachment, also resulting in a low structural reliability. 

For these reasons, an appropriate repair methodology for restor- 


ing the functionality and recovering the maximum reliability has to be 


identified. The joint reliability of the structure—the repair segment, 
the attachment, and the parent structure— must be simultaneously 
considered by computing the reliability of each part associated with its 
respective stress states. The joint reliability of the structure can be 
used as a parameter which characterizes the overall effectiveness of 
the repair. Based on this parameter, any general configuration of 
structure and repair can now be evaluated and optimized. 

This investigation identifies and outlines the principa! steps asso- 


ciated with such an analytical optimization of general structural repair. 


II. STRUCTURAL REPAIR APPROACH 


A. GENERAL CONCEPT 

A vital attribute of composite materials is structural redundancy. 
In order to identify the appropriate repair methodology for composite 
structures. the underlying redundancy musi De considered: 

The objective of repair is to restore the redundancy; even ai 
expense of increasing the stresses. In other words, our goal is to 
restore the functionality and to maximize the recoverable reliability of 
the whole structure. The effect of the repair on the entire Stric MEKi 
needs to be considered (i.e., the redistribution of the stresses eam ۱ 
by the damage, and later the repair itself which may transfer the fail- 


ure to a different location). 


B. FUNCTIONALITY AND RELIABILITY 
To restore the functionality, elements such as the type of Air cman 
component and aerodynamic (smoothness) need to be examined. 
From the strength theory for composite materials, the main fac- 
tors that are peculiar to laminate strength analysis (neglecting the 
thermal effects) are: 
° laminae strength 
* laminae stiffness 
* laminae orientation 
* laminae thickness 


e stacking sequence 


Thus, for restoring the reliability (redundancy), we need to opti- 
mize the patch stiffness through the materials selection and the lay-up 
orientation (including stacking sequence), together with different lap- 
joint configurations (scarf, step, or external patch). Moreover, the 
optimization objective is to achieve equal reliability in the parent 
πο ατα, the joint (the adhesive), and the repair patch, thereby 
achieving the highest possible total reliability for the entire structure. 
Upon the development of analytical computation methods for the joint 
reliability of the entire structure, the optimal repair configuration may 
be analytically sought by "iterating the choice of material (including 
Ius the orientations of the repair patch, and the joint 
configurations." Figure 2.1 is an example. 

A composite material structure with a repaired area is subjected 
to a complex loading, as shown in Figure 2.1. To evaluate the reliabil- 
wm ihe Stresses in each Spatial location need to be determined 
through the stiffness and material iterations. 

This example is a boundary value problem, and for that the linear 
elastic field equations must be satisfied (see Ref. 2). 

Due to the external load (Ti-surface traction vector), the stress 


boundary equation is considered, 


ΠΠ μιν Ζ 291) 


where (Sij) is the normal and shear stress components and (vj) is the 


unit outer normal vector to the surface. 
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Figure 2.1. Typical Repaired Structure 





(c) 


Due to the kinematics, the strain-displacement relations are 
defined (for small displacements): 


نے 
προ (2.2)‏ )911 * ,0,0( = ;€ 


where (ei) is the strain components and uj is the displacement. 


Substituting the strains into the constitutive relation, 


O =c E Tm 
νη (۰2-2. )[ 


where ciim is the stiffness matrix. 


Memce, the resulting Stresses. are determüned by satisfying the 


equilibrium equation, 


0,0, t X,7 O i, j = x, y. Z (2.4) 


where X; is the body force (usually ignored in the analysis). 

Now, recall from the constitutive relation (equation 2.3) that the 
stiffness matrix (Cii) characterizes the material property. In contrast 
to isotropic material. where the stiffness consists of two independent 
constants, and more importantly, the stresses are continuous. in an 
orthotropic composite material the stiffness matrix consists of nine 
independent constants, which reduce to five for special orthotropic 
orientation (because with respect to material properties Z and Y axes 
are equivalent). In a laminate, when the material principal axes are off 
the structural axes and there is symmetry about Zz - O, the stiffness 


matrix consists of 13 independent constants. Moreover, the stresses 


are continuous across the thickness (t), and the stiffness is a function 
of the lamina orientation and the stacking sequence. 

For a thin lamina under plane stress, the constitutive relation can 
be expressed as [Rel. 2]: 


The stress-strain relation is 


Ox Qi; Qi 9316 ہج‎ 
Cy کے‎ Qi; (323 (326 € y 
9 xy (316 Qe 66 E xy (2.5) 


where the Qg theis ٣٥٥ stiffnesses, are the Cijim stiff- 
nesses specialized for plane stress condition in an orthotropic mate- 
rial, and written in contracted notation (i, j = 1, 2, 6). 

The stress-strain relation can be related to the kí? lamina, and 


later (for FEM) to the nth element, 
Wake, = [Qs] 1 
i = μυ ο (2.6) 


From the laminate plate theory [Ref. 2]. the resultant externa! 
force acting on the laminate shown in Figure 2.2 is obtained by inte- 
gration of the stress in each lamina through the laminate thickness. 


For example, 


ና 
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where Ny is the force per unit length of the cross-section of the lami- 


nate. Similarly, My is a moment per unit length. 





Figure 2.2. In-Plane Forces and Moments 


Finally, by substituting equations 2.5 and 2.2 into equation 2.7, 


and summing the forces (N;;) and the moments (Μη) for the whole 


structure, the following relations are derived [Ref. 2]. 


λα ٦ 
ΠΕ ο ο κ, (2.8) 


where 
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εἰ, and ki are middle surface values (not functions of z), €; is the 
strain component, kj; is the curvature, and wo is the displacement of 
the middle surface in the z direction. t is the thickness and zy is 
defined in Figure 2.1(c). In equations 2.8 and 2.9. Ἆς- 15 ۶۷۷٦ 
extensional stiffness. Brs is called coupling stiffness (Brs = O if the 
laminate is symmetrical), and Drs is called the bending stiffness. 

Therefore, the stiffness in composite material is not depende mi 
onlv on the material stiffness (Qrs) but also on the laminae thickness 
(t = Zk — Zk - 1) and the stacking sequence (zy). The stiffness matrix 
(Qrs) changes due to the orientation of the laminae as it corresponds 
to the structure coordinates. 

In short, to iterate the material stiffness for determining stresses, 
numerous independent parameters are involved. The parent structure 
configuration is given (the stiffnesses are known). But in addition, in 
composite repair. two configurations are concurrently being consid- 
ered— the parent structure and the repair (shown in Figure 2.1). For 
each of those configurations, the stresses must be evaluated and by the 
constitutive relation (equation 2.3) they must match in their bound- 
aries. Before repairing the damage, the parent structure with a "hole" 


configuration needs to be iterated, and then the repaired +٤٦ 
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configuration. Through the resulting spatial non-homogeneous 
stresses, it is impractical to a priori determine which configuration is 
optimal (i.e., as can be done for isotropic material). Thus, by calculat- 
ing the total reliability from the state of stresses, we unify all the iden- 
tified parameters into one single parameter for optimizing the repair 


configuration. 


C. CONSIDERATION 

In this investigation, we assume that the functionality of a 
repaired structure must be the same as before damage, so the repair 
effectiveness is measured by the post-repair reliability. To evaluate the 
reliability, comparisons of the structure before damage, before repair, 


and after the repair need to be done (see Figure 2.3). 


Parent Damaged Repaired 
slructure structure structure 


R BOR R'<R"<R' 
0 ۶ 


Figure 2.3. Repair Configurations 
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The merits of the different elements of design oí the repa 86) 
be consolidated into one single unified parameter and measured by the 
structural reliability. 

From the general concept of reliability, the number of probable 
failures is associated with the magnitude of the stress level (the ran- 
dom variable). For a structure of given dimensions, the higher the stress 
level the higher the probability of failure. The quantitative probability of 
failure can be computed if the cumulative distribution of failure is 
known. For a structure under non-uniform stress (spatial heteroge- 
neous stress), the structure may be visualized as being divided into 
contiguous geometric elements within which the stresses are uniform. 
The reliability of each element corresponds to the respective stress 
magnitude of that element. Equivalentlv, the magnitude of the stress at 
each spatial location of the structure can be thought of as being 
mapped into a respective point on the stress space (by appropriate 
stress analvsis). If all the stress magnitudes are within the domain 
bounded by the strengths, then the entire structure is consider CHES 
be safe, as shown in Figure 2.4. 

The mapping of the physical spatial location to the stress domain 
is carried out in stress analysis. The definition of the strength in the 
stress domain is a failure criterion which can be defined either de- 


terministically or probabilistically. 


1 





Figure 2.4. Reliability Criterion 


1. A deterministic failure criterion can expressed in terms of the 


scalar stress components, which are themselves a function of the 
spatial location of the structure. 


f(0,(x,))«1 safe 
6 اج‎ fail (2.10) 


where o; is the stress component in tensor notation (i = 1, 2, 8) 
and xy is the spatial coordinates within the structure. F(.) is a 
function representing the contour separating all the safe stress 


|. 


states from all the non-safe states. For example, for coupled 
strength. F(.) may be expressed in terms of a polynomial crite- 
rion, as described in Reference 3. 


2. A probabilistic failure criterion expressed in terms of scalar 
stress components as random variables. The magnitudes of the 
scalar stress are themselves a function of the spatial location of 
the ٣٦ 


5 (2.11) 


where R is the reliability of a physical element located at spatial 

coordinates xk within the structure, under stress components 

(G;j)k, ek and Bk are the statistical parameters of the probabilistic 

failure model associated with the component stress oj. The 

entire structure is safe if the joint reliability for assemblage of all 

the spatial locations of the structure is greater than or equal to 

the reliability specification. 
D. TRADE OFF 

In the spirit of the consideration above, the joint reliability is 
consorted measure of the all the parameters which need to be opti- 
mized for a given repair. To explore the effects of the repair. s mm 
include the size of the damage area (after the clean-up) and the patch 
type (material, stiffness, and type of joint), trade-offs between stiffness 
and strength of the specific parent structure, the joining segment, and 
the repair patch need to pe ٣٣٣ 

If, when the repair design is completed, the structural reliability 
is not met, the design parameters need to be iterated to an achieve an 
acceptable solution. This iteration involves the redesign of the struc- 


ture/ repair configuration and the selection of new material and/or a 


new lay-up configuration. 
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III. STRUCTURAL RELIABILITY (CONCEPT) 


A. GENERAL CONCEPT 

In Chapter II, we developed the repair design concept as it 
relates to reliability. In this chapter, we will explore the specifics of 
stresses and strength calculations and examine their relationship to 
reliability. This relation allows us to compute reliability based upon the 
failure criterion which is most appropriate to the specific composite 
under consideration. 

Historically, formulation of many failure criteria was intuitively 
taken to be mechanistically based on failure mechanisms, but in fact 
the criteria are phenomenological (because frequently, the plane of 
failure does not necessarily correspond to applied stress). Phenome- 
nological failure criteria represent (describe) failure states but do not 
explain failure causes (or mechanism). There is no uniquelv "correct" 
formulation. and there is no a priori known formula. 

The failure modes observed for composite materials are quite 
unlike those observed for isotropic homogeneous materials because of 
their anisotropic and heterogeneous nature. Nevertheless, it is essen- 
tial for design to predict the capability of the structural element under 
a complex loading. Judgment must be made as to whether the theory 
(analytical results) adequately represents (for a wide class of problems) 


the physical phenomenon (i.e., experimental results). 
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In the design of repair. we are looking at the parent structure 
characteristic and the repair characteristic as a whole system's ability 
to withstand a given load. Then, by the reliability characterization we 
determine the stress level (coming from the applied load), for which 
there is an acceptable number of probable failures. 

In a composite material, the lamina is the basic building block 
which consists of fibers imbedded in a matrix. In composite stris 
especially in repair areas, the stresses are spatially heterogeneous, 
both in forms of stress components and in their magnitudes. More- 
over. because the lamina is composed of numerous fibers and 
possesses inherent redundancy, the composite failure proce- -Sii 
sequential! (in stress level). At the same time, the actual behavior of the 
loaded structure often exhibits considerable data scatter in measured 
strength due to processing variability (which is unavoidable in large 
composite structures). 

(Anisotropic) failure criteria are mathematical functions (analvti- 
cal. numerical. graphical) which define the stress levels below which a 
Structure is safe. If the measure of safety is absolute, the failure criteria 
are deterministic: otherwise. the failure criteria are probabilistic Me 
magnitude of stress (one component for 1-D; six components for 3-D) 
at each spatial location of the structure is mapped into 3 ۴٦ 
point in the stress space. If all the points are within the domain 
bounded by the strengths, then the entire structure is safe. Ideally 8 
failure criterion should capture the physics of the specific failure pues 


cess and the sequential events of failure. 
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Therefore, in order to calculate the reliability of a structure with 
non-uniform stress, the state of stress at every spatial location needs 
to be computed. The influence of these stresses on the respective 
reliability in those specific spatial locations can then be evaluated by an 


appropriate probability model. 


B. FAILURE MODEL—WEIBULL DISTRIBUTION 

In 1958(?), Rosen recognized that the mechanism of composite 
failure originated from the statistical variability of the fiber strength 
[Ref. 6]. The failure process was observed to be sequential, initiating 
from failure of weak fiber sites; the load originally carried by the weak 
fibers is transfered to the adjacent strong fibers, thereby postponing 
immediate failure. Upon additional load increase, the local micro- 
redundancy is insufficient to sustain the stress concentration, leading 
to catastrophic macroscopic failure of the composite. Harlow and 
Phoenix [Ref. 4] refined the probabilistic modeling of sequential failure 
events and arrived at the "chain-of-bundle" model. The composite 
structure is modeled as a chain of bundles in series. Each bundle has 
length 6. generally referred to as the "ineffective length." The ineffec- 
tive length is a function of the matrix stiffness and the interfacial 
bound strength; it is analytically estimated to be an order of magnitude 
greater than the fiber diameter (the diameter of graphite fiber is 
around 5 um). Each bundle is assumed to be independent and identi- 
cally distributed, and the cumulative distribution function (CDF) can be 
expressed in terms of the fiber statistics. Harlow and Phoenix have 


also derived, by analytical recursive formulation, the statistical 


! 


strength behavior of large composite from small bundles. Since the 
bundles are in series, the weakest link function is applied. For fila- 
ment strength, the Weibull model is appropriate for the characteriza- 
tion of the strength statistics. For the composite, the strength may be 
considered as piece-wise Weibull, where the scale and shape parame- 
ters are modified by the Harlow-Phoenix chain of bundle calculations. 

For this investigation, we consider each spatial element as 
independent and ignore failure sequence and subsequent load 
distribution. This is equivalent to stating that if one point in any 
direction fails, the whole structure fails (failure can occur in any element 
of the structure, in either the parent structure or in the repaired 
portion of the structure). The reliability under this condition can be 
calculated given the failure distribution and the failure criterion. 

We use the Weibull distribution with parameters adjusted by the 
"chain-of bundle" model. The two-parameter Weibull distribution is 


expressed as: 


For any ith spatial location (considered in our calculation as an ele- 
ment), the failure distribution function (CDF) is: 


Ε(σ (X) B) سے‎ 1 - expl- (5) } (3.1) 


where 


Sij— stress component at location x; 
p,— the location parameter, representing the strength under the 


respective stress component (oj). It is dependent on the size of 
the element. 
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phe shape parameter, representing dispersion of the failure 
function. It is also dependent on the size of the element. 


Anen, the reliability at each point, o (X) 


7۹۹۹۹۸۷ 9) (3.2) 


Cis expl- 631] (3.3) 


Assuming that the failure (Bj) event is independent of all spatial 


locations, the total structure reliability (RT) can be represented by: 


R_=R.R,---R 


T 12 K 


K 
Rr = IIR, 
1 


where (Rk) is the reliability at each spatial location. 


C. SIZE—EFFECT 

In our approach, for calculating the state of stress we are using a 
finite element method (FEM). For that we assume that a configuration 
m aired structure is partitioned into geometric meshes n of size 1 
πιο of size 2. The intrinsic strength of each element may be differ- 
sehe choice of element size is arbitrary, and because the same 
structure is under consideration, the total reliability must remain the 


SACE. 
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where 
RT- the total reliability 


ΕΙ, R2- The element reliability (assume equal for simplicity) of each 
meshing. 


n— number of element for size 1 
m- number of element for size 2 


By Weibull distribution, 


exp|- 5 d = exp|- 8.) | ο) 
1 2 
We make the simplifving assumption that the two different sizes 


are not drastically different so that the same shape parameters a can 


be used. It follows that: 


)3.7( وت 


1 


3/5 


Therefore, we can evaluate the reliability of each elementmas 
relating its length or area to the unit element. In our approach, the 
reliability is normalized by the smallest element. The smallest element 
is assumed to be the ineffective size. In short, we need to first perform 
stress analysis, which is in sections, then evaluate the reliability 


knowing the strength of the smallest element. 
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IV. REPAIR METHODOLOGY 


A. TYPICAL REPAIR REQUIREMENT 


The previous chapters discussed the rationale of adopting the 


joint reliability as a single indicator of the repair efficiency and the 


property of a probabilistic failure criterion which is required to calcu- 


late the reliability of each spatial element under a non-uniform 


complex state of stress. Furthermore, we identify the following 


requirements which are necessary to establish the appropriate repair 


methodology for composite structures with the goal of restoring the 


structural integrity: 


1 
2 


The entire structure needs to be considered. 


The patch stiffness through the materials selection and the lay-up 
orientation, together with different lap joint configuration, need 
to be optimized. The optimizing objective parameter is the joint 
reliabilitv of the structure, including the parent structure, the 
mi and the repair patch. The objective is not to maximize the 
reliability of the repair but rather to achieve highest possible total 
i liability for the entire structure. 


. In order to assess the necessity or the effectiveness of the repair. 


a single parameter is formulated. We select reliability as the 
parameter (refer to Chapter MI). 


. Using the reliability formulation, the parent structure reliability is 


compared to the damaged structure reliability before repair, and 
the repaired structure reliability. As mentioned previously, the 
ideal repair restores the reliability of the repaired structure. 


. It is known that conventional repair depends on the adhesive 


joint, which always has no redundancy or poor redundancy 
(smooth-joint versus step-joint). This weakness has to be taken 
into account. In our approach, the adhesive strength affects the 
stress distribution and limits the repair stiffness (see joint con- 
figuration in appendix). Thus, the joint strength is considered 


21 


with the stiffness optimization (lay-up design) of the repair seg- 
ment (i.e., weak adhesive limiting the transferring of stresses to 
the repair segment compliant patch) 

6. To simplify the calculation, and to reduce the number of parame- 
ters that have to be iterated in the same time, we assume the 
adhesive reliability to be equal to 1. In other words, we assume 
infinite joint reliability. 

7. The reliability of the structure is the product of the reliability of 
the individual elements (see equation 3.4). Moreover, to achieve 
the reliability of the element, it is required that the magnitude of 
the stress tensor not exceed the failure criterion (i.e., a spatial 
reliability function with stress tensor as the random variable (see 
Chapter II). 

B. CONFIGURATIONS FOR REPAIR 
Now, in order to utilize the appropriate procedure for designing a 
repair, we have to consider three separate configurations: 

l- The parent stricture 

2. The damaged structure 

ከ ርን 0 ን 
1. Parent Structure Configuration 

For the parent structure, a finite element model must be built 
to represent the original geometry, the lamination configuration, and 
the original material. The boundary conditions are determined by the 
actual external load and reactions (support). 

By knowing these data, the resulting stress tensor for each 
element can be calculated. Using the suitable failure criteria, the 
parent structure reliability (RO) can be evaluated. The parent structure 


reliability is the acceptable parameter reference to compare for 


optimizing the repair configuration. 
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2. Damage Structure Configuration 

We assume that all the facilities needed to evaluate the 
damage area are available (i.e., NDT). Thus, before removing the 
damage and preparing the area for the bonding repair, the clean-up 
geometry has to be optimized. Dependent on the size of damage, the 
spatial location, the repair capability, and the state of load in this area 
(considering the fibers' orientation), several configurations can be 
evaluated. Again, going through the same procedure introduced for the 
parent structure, the state of stress is calculated and thereby the 
reliability for each configuration, R1, R2, R3. Now the configuration 
with the largest reliability (say R3) is selected. 

3. Repair Configuration 

Starting with the R3 configuration, and assuming infinite 
joint reliabilitv, the joint strength is considered with the optimized 
stiffness. Thus, onlv the reliabilitv of the parent structure (outside of 
the repair zone) and the reliability of the repair structure (the patched 
zone) are considered for calculating the total repaired structure relia- 
bility. To optimize the repair configuration in the lay-up design, one 
must control the repair stiffness and redistribution of the stresses. 
Either one type of material can be chosen and only the lamination can 
vary (lamina orientation) or more than one material and lamination can 


be applied. 


C. FINITE ELEMENT APPLICATION 
It has been shown that to calculate the reliability, we need to ana- 


lyze the state of stress in all spatial locations— analytical boundary value 
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problems or finite element methods (FEM) can be used. Because of the 
complexity of the structure geometry, it is essential to use FEM (refer 
to ADINA code in Appendix A). 

The analysis of stress in a laminate of laminae is straightforward. 
For 2-D problems, we are using the laminate plate theory (plane 
stress). For 3-D problems, we are using generalized calculations for 
orthotropic materials. Assuming the laminae exhibits linear elastic 
behavior, the stresses seen to be a linear function of the applied loads. 
That is, if all laminate stresses in each laminae ({i.e., corresponding to 
an element in ADINA code) are known, then the stresses in each 
laminae can be compared with the lamina failure criterion (see Chap- 
ter II) and uniformly scaled to determine the load at which failure 
occurs and the spatial location. 

In accordance with the reliability and failure criteria, all failure 
criteria for laminates depend on the strength in the principal material 
direction, which likely does not coincide with principal stress direc- 
tions. Therefore, the strength in each lamina in a laminate must be 
assessed in a coordinate system that is likely different from those of 
its neighboring laminae. (This coordinate mismatch is one of the com- 
plications that characterizes the determination of the failure criteria.) 

Hence, for the interpretation of stresses resulting from the FEM 
calculation, all the stress tensors have to be transformed to the prin- 
cipal material direction per lamina (see Appendix C). 

Furthermore, in order to normalize the reliability and to make it 


geometrically independent, the area of each element has to be related 
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to an acceptable minimum area. The minimum area is derived from 


the minimum size effect. 
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V. METHODOLOGY FORMULATION 


The object of this chapter is to explore the effects of the damage, 
and finally the repair, on the stresses of the entire structure. Appro- 
priate simplifying assumptions have been represented in the the 
repair methodology (Chapter II). 

We wish to capture the essence of the underlying mechanics and 
physics of our problem. The field equations in solid mechanics are 
specialized for calculating the state of stress in the repaired structure. 
Then, the acceptable probabilistic failure criteria, expressed as a 
mathematical function which defines the stress level below which the 
structure is safe, is utilized. 

Therefore, in the context of preliminary sizing and design trade- 
offs of the repair. and also for confirmation of output of more accurate 
analysis carried out by the computer code. we represent the repair 
methodology by the following formulations. 

The formulations are classified by : 

l. The physical space- a line, a plane, a volume 


2. The stress space-— one component, three components ٦ 
components 


3. The spatial distribution of stress within the structure— homoge- 
neous or heterogeneous 


4. The definition of safety— for our case, probabilistic failure criteria 
Whereas in one-dimensional theory (1-D uniaxial stress), each 


scalar stress is mapped into a respective point on a line (1-D) in the 
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stress space, under planar combined stress, the magnitude of stress 
(three components) at each spatial location of the structure is mapped 
into a respective point on the 3-D stress space. Under general com- 
bined stress, the magnitude of stresses (six components) at each spa- 
tial location of the structure is mapped into a respective point on the 


6-D stress hyperspace. 
Stress interior to domain bounded by the failure surface — safe. 


Stress exterior to domain bounded by the failure surface — not safe. 


A. ONE-DIMENSIONAL EXAMPLE 
In order to understand the methodology principle, a simple stress 
analysis, consisting of a scalar calculation for single component stress, 


is shown in the following example. 


A long segment (shown in Figure 5.1) of non-uniform cross-sectional 
area (A(x)) is subjected to a uniaxial tension load (P). Then. to cal- 
emirate the reliability: 


1. The long segment is divided to small and unequal elements. 
because of the usual complexity of the repair geometry (geometry 
mapping) and because of different stress magnitude distributors 
(i.e., the segment is spatial stress heterogeneous). 


2. The stresses (o(x)) are uniaxial stress (1-D) and must be evalu- 
ated in each spatial location. Thus, because of uniaxial stress (one 
component, 1-D), the stress analysis is simple and represented 


o(x) = _P 
A(x) (5-1) 










1 


Figure 5.1. Long Segment, Non-Uniform Stresses 


Therefore, because of the series process assumed for composite 
material structure (refer to Chapter II), the total reliabilitv of the 
structure (RT) is presented as a function of the stress and the 
2۱611161717572 ٦ 


RR IRO (5.2) 


where 
n— the number of elements 
Ax 21 117 elen ٣ 


Furthermore, using the Weibull distribution function (weakest 
link of chain), the total reliability is expressed as, 


es (5). sa 


. To compare the structure configuration )1.6۰,۰ 1707 nB 
different type of repair), and to avoid size effects on the strength 
of the material, the reliability must be geometry independent. 
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Thus, the reliability is normalized by the smallest length, as 
follows. 


Using the weakest link rule [Ref. 4], and assuming the shape 


function is constant (in composite material a@ - 5-10), the size 
effect is represented by, 


ፔን) (5.4) 
B= aay (5.5) 


ποτε 
0— the length of the smallest element 
p— the mean strength value of the smallest element (evaluated 
by experiments) 


p;— the mean strength value of the ith element. 


Then, by substituting equation 5.5 into equation 5.3, the total 
reliabilitv is obtained as: 


[E Hexp - (S (5-6) 


. Finallv. a simulation needs to be performed to demonstrate the 
probability of different loading configurations. This simulation was 
carried out via a statistical analysis using the spreadsheet Micro- 
soft Excel on a Macintosh computer (see Reference 9). There, 
the strength is evaluated by the Weibull distribution function 
(weakest link), and if any stress is larger than the suitable 
strength, then the structure fails (1— fail, O— no fail). 


Thus, after several runs of the simulation, we can observe that the 
location of highest stress has the highest probability of failure but 
not a certainty of failure. In other words, the failure is not con- 
stant due to the spatial location— it is changing according to the 
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probabilistic failure criteria. The resulting reliability distribution 
is compared so the configuration of the structure loading can be 
optimized. 
B. TWO-DIMENSIONAL EXAMPLE 
The first approach that corresponds to the macro-mechanical 
behavior of a composite laminate is the 2- D stress calculation using 
the classical lamination theory [Ref. 2]. Thereby, under the same prin- 
ciple as 1-D calculations, 2-D calculations are performed, considering 


three stress components: 


πο σα, 


A typical analysis of a thin plate with an elliptical hole subjected to 
tensile load in a non-principal direction is described. This analysis 
models the damaged area of a structure after "clean-up" and before 
the bonding process of the "patch" repair. 

Again, going through the same procedure as the 1-D example in 


order to calculate the total reliability, 


l. The physical space considered is the space area of the entire 
plate (laminate). The space is divided into unequal areas. The 
rectangular rather than the polar division is preferable because of 
the orthogonality of the composite material, and because it is 
easier to compare the structure reliability as a spreadsheet. 


2. Assuming an unequal area for each element, the stress analysis is 
performed. 


Note that for using reliability equations, the stresses in each spa- 
tial location must transform to the principal material direction, 
as explained in Chapter II. For this, stress transformation equa- 
tions are used [Ref. 2]. 


Now, having the stresses in each spatial location, the same pro- 
cedure as in the 1-D example is utilized. We must remember that, 
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here, the weakest link rule relers to each component of stress in 
each element (three components). Therefore, the weakest link 
ance ers to each clement in the entire structure (plate). 


3. To normalize the reliability, the lengths in the previous example 
(equation 5.5), are replaced by areas: 


1 


e O 


k.i 





k— relates to the element notation. 
i— relates to the stress component (iz1,2,6). 
Ar — the smallest element area. 


Ay ;— the kt! element area. 


4. Therefore, because of the weakest link approach, there are two 
options to calculate the total reliability (Rr). 


Option 1: Calculate the reliabilitv of the whole structure with 
respect to each component in each element (shear or normal 
stress), then calculate the total reliability, as shown below. 


The reliabilitv of the k!? element due to the it? stress compo- 
nent is 





x 
R.,= exp -[ | i= 1, 2, 6 (5.8) 
ο 


Note that the stresses are in contract notation. 


The whole n elements reliability with respect to the ith stress 
component! is: 


R.,=LIR,, LL (5.9) 
k=1 


Finallv. the total reliability is: 
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R4 - ]L ÍR,, i= 1, 2.6 (5.10) 





[NE T ITIR,, - 1111 - የብ] 


i=] k=] i=] k=] Du 


Option 2: calculate the reliability of each element and then cal- 
culate the total reliability as follows. 


The reliability per element (Rg.x) is the product of the stress 
component reliabilities (Ry ;j— expressed in equation 5.8). 


Thus, 


6 
ር Πε, Γι ο (5.12) 
i=] 


Then, the total reliability is, 
Με - (5.15) 


And, after substituting equations 5.10 and 5.14 into equation 
5.15, the total reliability is: 


n 6 a ። 
ርን η i= 1577 
8, -1111፡58-(፤ | Κι ክብ (5.14) 
k.i 


ዜ=] 1=1 


Although the total reliability must be the same for any option, a 
difference exists between the two options and is obtained by the 
interpretation of optimizing the repair. Through Option 1, we can 


visualize the effect of a specific stress in a specific direction (i.e., 
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improve fiber, if the total reliability— [Rm ;— fiber direction], is small). 


Or, through Option 2, the weakest element which has the lowest reli- 
ability (R£ k) can be investigated. 


C. ANALYSIS IN THREE DIMENSIONS 

The next obvious approach is to extend the analysis to a laminate 
in 3-D physical space considering combined and spatial heterogeneous 
stresses (six components). Obviously, the stress analysis is more 
complex and cannot be solved "by hand,"—a numerical analysis 
method must be used (see application of ADINA code in Appendix A). 
Moreover, following the same procedure that is used in the 1-D exam- 
ple and in the 2-D example, the total reliability is calculated. 


1. The state of stresses resulting from the 3-D analysis (generalized 
orthotropic laminate— Reference 2), must relate to the mean 
strength known for laminae (from experiments). In other words, 
the stress must be calculated per element in the lamina. There- 
fore, unit volume will have an identical thickness for each ele- 
ment Hence, the area normalization in 2-D (see equation 577) is 
sufficient also for the 3-D calculation. 


2. Furthermore. to calculate the reliability, instead of three compo- 
nents in the 2-D example, we consider six stress components: 


19,.0,,0,,0,,0,,0,1 


For the first option mentioned in the 2-D example, equation 5.10 
extends to: 





ΘΠ 
R, = exp (2) | ο 8 (5.15) 
κ, 
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The whole n element reliability with respect to the ith stress 
remains the same (see equation 5.9). Hence, the total relia. 
IS: 





اتا 
መ‏ 


.-1[[ኮ.. -1111፡፡8- e k= 1 


i=] k=1 i=l k=] 


For the second option mentioned in the 2-D example, equation 
5.12 extends to a product of six instead of three, as shown below. 


Rp EN MMC (5.17) 
Thus, similarly to the 2-D example, the total reliability is: 


š i21,2...6 
s. -Tllles - mm 11 k-12...n δι 
Kat 


πιει 





SUMMARY 


. The 1-D example is only a simple presentation of the methodol- 
ogy concept used for calculating the total reliability. 


. In the 2-D example, the whole structure analysis is simplified 
using laminate plate theory. This approach is used in practice to 
do preliminary calculation and sizing. and to check the results of 
2-D and especially 3-D numerical analysis carried out by com- 
puter codes (see Appendices A and B). 


. For the 3-D example, the state of stress in each spatial location, 
considering the generalized orthotropic laminate, is numerically 
calculated using FEM by computer a code (i.e., ADINA). 


. Comparing 2-D to 3-D as the generalized orthotropic laminate 
analysis, in 2-D analysis the interlaminar stresses are ignored. 
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VI. CONCLUSIONS 


This work has proposed a primary approach for the designing of 
composite material repair, using reliability as a single parameter to 
indicate the most effective repair. Finite element analysis was used to 
calculate the state of stress at any spatial location (non-uniform and 
combined), which in turn was used to calculate the resultant reliabil- 
ity. This methodology can be used to assist maintenance personnel in 
the selection of the optimal repair configuration. 

The Weibull distribution function with modified shape and scale 
parameters was used for estimating the reliability, based on the obser- 
vation that composite material failure is adequately represented by the 
weakest link chain-of-bundle model. The best selection of failure cri- 
Fri πας to be confirmed by experiment. 

Simulations for different patches under a specific loading or for 
Specific patches under different loading need to be included in the 
future extension of the repair optimization method. Moreover, an ade- 
quate number of samplings (experiments) needs to be performed in 
order to estimate the underlying Weibull parameters and then to 
approve the appropriate failure criteria. 

Only bonded repairs with infinitely strong adhesives have been 
considered. Further studies may be extended to apply this approach to 


consider adhesives with finite reliability and bolted repairs, which are 
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considered as 3-D reliability enhanced (small transverse fibers are 


included under this category), and to scarf joint geometry approval. 
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νι σοι ۹ 
DESCRIPTION ΟΕ ADINA 


A. GENERAL DESCRIPTION 

ADINA is a finite element code for automatic dynamic incremental 
nonlinear analysis. The ADINA system is composed of four main pro- 
grams: ADINA-IN, ADINA, ADINA-PLOT, and ADINAT. Version 2 (Jan. 
1989) was used for this thesis. The sequence of execution of these 
programs is shown in Figure A-1 [Ref. 6:pp. 1-2]. 

mon Our investigation, the ADINA code was used to conduct the 
stress/strain analysis (3-D/2-D) of an advance composite material 
structure. Thermal effects were ignored, so the ADINAT program (for 
heat-transfer analysis) was not used. 

The ADINA-IN program generates the input data for the ADINA 
program and has the ability to graphically check the model. All or part 
of the model can be drawn from any view desired, including element 
and nodal numbering. 

Coordinates can be input in a local or global coordinate system 
using cartesian, cvlindrical, or spherical coordinates. 

ADINA is used for displacement and stress/strain execution, pro- 
duces the printout file for displacement and stress/strain, and also 


٣٠۰٠٠٠٢٢٣ the porthole file for ADINA-PLOT. 


ST 


COMMANOS 
OATA LINES 


0 πελο ειτε 


ADIHA - IN 
OATABASE 


AO!NA - IN 


AOINAT INPUT 
AQINA INPUT 





- 
------------ --- — no 


AOINAT TEMPERATURES 


ADINA 


RESTART 
INFORMATION 





PORT HOLE 


I 


AQINA-PLCT ۸0۱۷۸1۲ 
OATABASE 





Figure A-1. Sequence of Execution of ADINA Programs 
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ADINA-PLOT performs the post-processing of the output data 
from ADINA. It also has graphics capability, allowing a view of both the 
deformed and undeformed shapes, and can do vector plotting of 


stresses. 


B. COMPOSITE MATERIAL IMPLEMENTATION 

The program defines the parameters (modulus of elasticity (Ej). 
Poisson ratio (vy) for the orthotropic linear elastic material model 
[Ref. 6:pp. 5.5-6). The material defined is assigned an identifying 
ci which can be referenced by several or all of the element 
٣١٠٠٠٠٢ ie., parent structure elements, repair patch elements). The 
material directions correspond to the global coordinates of the finite 
ament model [Ref. 6:pp. 5.5-7). 

In ADINA, the composite material model can be used by several 
٠۰۰٠۱٠٣ ٦٦065: 2-100 solid element, 3-D solid element, plate element, 
ell element [Ref. 6:p. 5.5-3). In this thesis, the 2-D element and 


3-D element type were used. 


C. VERIFICATION 

Part of the establishment of the repair guideline is the verification 
sme finite element code. This is done by modeling a problem for 
which the results are known and comparing numerical with these 
known results. 

The verification includes the geometry verification, material axes 


definition, and their transformation to the global coordinates, and 
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compares the 2-D (plane stress) solution and the 3-D solution with 
theoretical solution done by hand. 

A simply supported anisotropic square plate, under uniform ten- 
sile load along one edge, is considered as a verification problem and 


discussed in Appendix B. 
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ھ٤‎ 
VERIFICATION PROBLEM 


A. PHYSICAL PROBLEM 

A simply supported anisotropic square plate, under uniform ten- 
sile load along one edge, is considered. The plate is built of eight plies 
and is symmetrical, laminated once with stacking sequence [902/759]. 


and once with stacking sequence [904]s. 


B. OBJECTIVE 
To verify the performance of the 2-D solid plane stress element, 
and then the performance of the 3-D solid element, when the material 


is anisotropic. 


C. FINITE ELEMENT MODEL 

For 2-D analvsis, eight-ply thickness plate (laminate) is consid- 
ered. The mesh consists of three 2-D solid elements (nine nodes). An 
orthotropic linear elastic material model is used with the principal 
material axes a and b, making angle beta with the global coordinate 
axes X and Y, respectively. The laminate stiffness is calculated bv the 
mus bi C program [Ref. 7:ch. 7] from the specific sequence and 
orientation of the laminae in the laminate. 

For 3-D analysis, each ply (lamina) is considered separately across 


the thickness. The mesh consists of 12 3-D solid elements (27 nodes). 
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The same orthotropic linear elastic material as in the 2-D model is 
used. The lamina stiffness is used directly in the finite element model. 

For any model, through ADINA execution, a printout file is cre- 
ated. The printout file generates a listing of the displacements and the 
stress/strain for each element. Only the stress/strain results of nodes 
number 9, for the 2-D problem, or node number 21, for the 3-D 
problem, are printed. Node 9 and node 21 are located at the element 
geometric centroid and represent the element mean values for calcu- 


lating the reliability. 


D. SOLUTION RESULTS 

The solution results are as expected. For the lay-up of [904]s, the 
stresses in the applied loads direction are the same as the applied 
loads. No interlaminar stresses are shown (tend to zero). The results 
of 2-D analysis are in good agreement with the results of 3-D analysis. 

For the lay-up of [902, 759], in 3-D analysis, the layers in 90° ori- 
entation are carrying more stresses than the layers in 75° orientation, 
and the interlaminar stresses are shown (they cannot be neglected). In 
2-D analysis, because of plane stress assumption, the interlaminar 
stresses are not shown, and the stresses have the same value as 


before— they are representing the mean stress. 
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P Bc ch Be Ας e ہار ہاب بر ہل ہاو پار چاو بار‎ he hh Ας Ας e e hk hh hh hh hh he hh Ak 
PROGRAM TRANS2 

Be Ae i hh kh he hh e hh Ας e e hh e پل‎ e x e hh he he hh κά κά κ e Ak 
INTEGER NLA,NEL 
PARAMETER (NLA=6,PI = 3.1415927 ) 
PARAMETER (ROWS=100,COLS=6) 
REBE~G A(ROWS,COLS),B(ROWS,COLS) 
Ree Mol ,;MOZ,NULZ,NU21,EX,EY,EXY 
REAL THETA(NLA) ,M,N 
٣۰۰۰٠٠۰ Ol1,912,022,066,G12,E1(NLA) ,EZ(NLA) ,E6(NLA) 
CHARACTER^40 ,FNOUT 


EE (5,*) ' OUTPUT FILE,NAME: ' 

READ (5,'(A)') FNOUT 

ΙΙ. “INPUT MODULUS & POISSON MOl,MO2,NU21,C12 ?' 
۲٢١۰٠۰٠٢٠٢٠٢ ٦ , ۸(۰ 101,102,102 72 


NU12 = M02*NU21 / MO1 
Il =(MOQ1 / (1=NU12*NU21) 
Q12 = Q11*NU12 
ረረ 9 1-NU1LZ*NU2)) 
Q66 - C12 


Pere (5,%) "INPUT STRAIN EX,EY,EXY ?' 
Ema 5^) EX,EY,EXY 

OPEN (UNIT-13,FILE-FNOUT,STATUS-'NEW') 
MITE (13,x) “INPUT : ' 

NIE 5, 10) ΕΕ ΕΝΥ 

ο τι (13,10) EX,EY,EXY 


11 5 x) “INPUT THE ORIENT. OF EA. ELEM. IN DEG. 7?’ 
Pewee (13,*) “INPUT: THE ORIENT. OF EA. ELEM. IN DEG. 
DO 50 I=1,NLA 

ΠΕ (5,20) I 

BEAD (5,*) THEIA(I) 

۶۰٠۰٠۰۶٢٣٦3 ,25( ۲1197513 11[( 

οσο TIETA( )#*PI 7 180,0 
50 CONTINUE 


. 
r 


* 

EINE (5,*)'OUIPUT :' 

Meme (13,*) '"OUIPUT : ',FHOUT 

۰۰۰.۰٠٠ 60٣۳ "NLA' ,' 51 52 53,56’ 55’ S4 

BEIPO(15,60) 'HDLA','81','S2','$S3','S6','S5','S4' 
* 

NEL-0 
DO 90 J=1,3 
DO 80 I=1,NLA 
NEL=NEL+1 

* STRAINS TRANSFORMATION 


M=COS(THETA(I)) 
ΕΙΝ ΤΠΕ ΑΕΙ) 


F b i) = حم‎ ተየ ጠለ B EYAI*A2 + EXYAMAN 
|| .=... ተ እጋ بې‎ EYAMAK2 + EXYAMAN 
ο ወ. μμ + EYA2Z*MAN + EXVYACMAA2D >= N&*2) 


ኣ SIRESSES CALCULATION 
ΙΙ eel tel Gl) + O12Z4E2Z( 1) 


A(NEL,2) = Q12*E1(I) + Q224E2(1) 
A(NEL,6) = Q66*E6(I) 

A(NEL,4) = 0. 

A(NEL,5) = 0. 

A(NEL,3) = 0. 
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80 CONTINUE 
2n CONTINUE 
DO 100 T= NEL 
WRITE (13,70) NEL,A(1,1),A(1,2) Αι ο ο. NIIS 
100 CONTINUE 
ዴ 


CALL BETA(A,B,NEL) 
WRITE (5,%) “INBUT ALFA 
READ (5,%*) ALF 

WRITE (13,40) ALF 

CALL ALFA(A,B,NEL,ALF) 
CALL RELI(A, B, NEL) 

CALL ELRE(B,NEL) 


0 FORMAT (//,1X,'STRAINX= ' ,۲11.3, 87, 6 ات لگ د‎ 
& 9X,'STRAINXYs ',E11.3) 

20 FORMAT τη ΤΡ τ; 0) 2) 

25 FORMAT (IX, TELA TIS STS ToD 


40 FORMAT (//,]X,'SHAP FUNCTION (ALFA) 5 Ες 1 ۴1 6 
60 FORMAT ο 110,6 ۸ ۹ κ οι OI 
70 FORMAT ረ To 0۹ 
x 
STOP 
END 


Ck o e o fe hh AA ch A o II A 
SUBROUTINE BETA(A,B,NEL) 

BoA AR AR AR A RCA OR hc ch A AKA 
PARAMETER (ROWS=100,COLS=6) 
REAL^8 A(ROWS,COLS),B(ROWS,COLS) 
REAL X (COLS ΓΤ; 


71۱+ "8 ۰۰7۷ ۰م 


WRITE (13210) “NEL, SRI 312 ον ο ο. 
DO 90 151 ΝΕ. 


X(1) = 15.0E4 

X(2) = 0.4E4 

፳(3) = 2) 

X(6) = 0.68E4 

X(4) = X(6) 

X(5) = X(6) 

πιο 5 IS ٤ 

۸648 ھ۔ نے 20ے 

تک الا 
LF τμ;‏ 
IF (A(1, 2) 61 0155 02) = 225‏ 
oe‏ لاه IF (ACI 3) DP O)‏ 


DO 80 J=1,COLS 
BIJ) = ۰  ) +۶ ۵0 

BO CONTINUE 

WRITE 113,60) Ι,ΒΙ1,1),ΒΙΣ,2) , 8٣" 07٦ 
90 CONTINUE 
10 FORMAT (/74,1X,Al10,G6 CA 1 4X) ٢ ل٤‎ OD 
60 FORMAT (7, 4X 6,6 ٦٣( 

RETURN 

END 


KkkkkAKkKkkhkkkkAkAkAKkkAKkKkkkAKkAkAkAkkAkA KAA 


SUBROUTINE ALFA(A,B,NEL,ALF) 
پا او يهر‎ coh e eh o o o f i oh e oh e he پر لو پلر لو لو‎ 
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PARAMETER (ROWS=100,COLS=6) 
REAL*8 A(ROWS,COLS) ,B( ROWS ,COLS) 


E SIRESS = SIRENGTH,RATIO RISED BY ALFA CALCULATION 


INE ΕΠΕΙ͂ VI V2, UI V6 νο, νὰ” 
ΓΟ 90 I=1,NEL 
DO 80 J=1,COLS 
AUI,J) - B(I,J)^*ALF 
BO CONTINUE 
DENIM 13,60) L7 AVE71),ACI2)7ECT, 3), ACI, 1 ,IINLFEZS) ; ACT, 4) 
90 CONTINUE 
10 መክ |፡//.11,አ130,6(።እ/,5ሯ) / 1]2፲,80('-"ነ) 
60 BHORMAT-(/,4X,I6,6E11.3) 
RETURN 
END 


e Bh K e dh e e I e ee b e e وياو لو باو يريو باو ياو باو‎ 
SUBROUTINE RELI(A,B,NEL) 

Te Pe Bebe K e Pe k k ee Pe e Be e e e e e e e he κά 
PARAMETER (ROWS=100,COLS=6) 
REAL^8 A(ROWS,COLS),B(RORWS,COLS) 


EE UN SIERSS TENSOR RELIABILITY CALCULATION (WEIBUL) 


NNEESCUIS.10) 'NEL',' R1','R2','R3','R6','R5','R4"' 
DO 90 I=1,NEL 
DO 80 J=1,COLS 
EIU)  EXPCSA(CI,J)) 
BO CONTINUE 
۲۰-۱ CEU), I;BCI,1),B(CI,2) ;BUC1, 2) , BCI,6), B(I, S), BCI,4) 
90 CONTINUE 
10 fete (//,142,À10,6(A7,472) / 12,580('-፦')) 
60 EMT (/,4X2,16,6E11.3) 
RETURN 
END 
& de Pe Pe # K # e he he de ja de jo he لو لو مو‎ pe he he he he de de باو بل لو لو‎ 
SUBROUTINE ELRE(E AWEL) 
Αν h ἐς خو بو باو بو بو باو باو باو لو باو‎ pe خو بو‎ de po he خر بل خو خو لو بل لو باو خو خو‎ 
PARAMETER (RO0WS-:100,COLS726) 
fees Bi rUwe ,COLS) ,ELe(ROnS, ,STH 





κ 
LS SUS RELIABILITI ChLCUL. &s STRUCTURE RELIABILITY CALCUL. 


DOE (12,10) 'MEL','ELM.RE. 
pen = 1.0 
DO 96 I-:1,5EL 
BERI) = 1.0 
DO 66 J=1,COLS 
EDD) = ELI) N ELI) 


80 CONC TMUE 
N'TIE (5,50) I ELR(I) 
NETPIE (15,60) I,ELKR(I) 
ome = SIR + ELFI) 
90 CONTINUE 
NELLE (5,70) STR 
RITE (13,70) STR 
10 ΤΠ 2510 / 1Z,25(”-')) 
60 ΙΙ Ιι'!,4;,16,Ε]1].3) 
70 ΙΙ, ΤΕ, STRUCTURE RELIABILITY = ',FE11,3 /1ች2,40('=')) 


APPENDIX G 
POST-PROCESSOR PROGRAMS 


We have the state of stress in the spatial locations defined from 
the FEM (ADINA code). The next steps for evaluating the reliability 
are: 

1. To normalize the size of the element. 

2. To perform, per element, the transformation of stresses, read 
from the FEM output file, to the lamina (composite material) 
principal directions. 

3. To compute the reliabilty based on a specific reliability model. 
This process is done by two programs— CLAREA and TRANSS or 

TRAN DS 

CLAREA is the program which reads in the coordinates as well as 
the element area defining the nodes. It then computes the element 
area and the smallest element area. 

TRANS3 is the program used for 3-D analysis. lt reads gm ٢ ا‎ 
stress components of each element as well as the element area and 
the angle theta, coresponding to the angle between the global coordi- 
nate system and the layer pricipal coordinate system. It then performs 
the transformation of stress to the principal direction. In addition, 
this program calls the subroutine BETA, which in turn calls subroutine 
ALFA, which in turn calls subroutine RELI, and which finally in turn 
calls subroutine ELRE. 
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TRANS2 is the program used for 2-D analysis. It reads in the 
strain components as well as the modulus of elasticity utilized in the 
repair or the parent structure, the angle theta, and the element area. 
It then calculates the stress tensor using the constitutive relation. 
From this point onward, the program is the same as the TRANSJ3 
program. 

These programs are illustrated by solving the verification exam- 


ple, as shown below. 


47 


2-D PROBLEM 
(Plane Stress 


[30415] 
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OE TT neur ؟‎ 


ት ታት ታታ 


s ነ! ANALYSIS OF UNISOTROPIC PLATE 


Kk 

ENISSUNITS LIST=8 L0G-7 ECHO-7 
CONTROL PLOTUNIT=CM 

COLORS BCODE=BLUE 

κ 

DATABASE CREATE 

WORKSTATION DEVICE=0 

* 


ES X IHV2D.IN ANALYSIS OF UNISOTROPIC PLATE ( 90,90 DEG )"' 
* 
ERSTER IDOF2000111 


ντου VOLUME-MAXIMUM P=STRAINS 
ΙΡ FORMATTED-YES FILE-60 
x 


COORDINATES 

ሽክ [FS NODE Y ὦ 
1 20. 20. 
DE 20. 20. 
3 20; ως 
EE 20. S 
5 ΖΩ z 95 
E IG. E 
1 20. 2207 
fee 20s 203 
9 0. m 

* 

* 

fee 123 / 9 

meme ls / 78 

x 


u I] ORIHOTROPIC EA=181.E05 EB=10.3E05 EC=10.3E05, 


NUABSOS0ISS.NUACSO.0159 NUBG=0.26 CAB=7. 27505 
۸ 


PLOUR TWODSOLID STRESS2 TABLES 
ΙΕ 2 4 3 EL1-1 EL2-1 NODES-9 
|” ስ ος 5 6 8 7 EL1:1 EL2-1 NODES-9 
STABLE 1 9 
72 
AMTES- EL THICK EE TABLE 

1 Creal 90 1 

2 On 90 1 


* 
LOADS ELEMENT 

1 2 -3.0E06 
ል 
EGROUP 2 TWODSOLID STRESS2 TABLES 
GSURFACE 3 4 6 5 ELl-1 EL2-1 NODES-9 
STRESSTABLE 1 9 


EDATA 

RIES EL THICK BEI TABLE 
1 Oni 90 1 

κ 

A 


VIEW ID-1 ZVIEH-1 YVIEW-0 ZVIEW=0 
EGZONE HAME-TWODSOLID 
1 
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FRAME THEADINGSU pn ٥٣ YSF=0.5 APMAX=22.5 /1/ን 
MESH ZONE-TWODSOLID NODES-11 ELEMENT-2 BCODE-ALL VIEW-1 
x 


ADINA 
۸ 


END 


0 


ረ. 1 ኮኮ: 5 ዐበ- εὐ 0 £O SALES ٣ 0 U> 


00 ’ 06- ed psp) 50+3000 00 ει οί ο- S0rd000t'0 ۰38٤٣ 0000 0 6 
ር 
ىےۓ‎ 38086 0> ٦0-1 0 ۲۲۳ O2 CO IIE UO 
00 ° 06- 6۱۲935606 O02 50+3000 0 6-0 0- 0 1۱4711) OZ 00+00 6 
I 
(SHLWOIGNI OS (SdI) 9۷134 111911013 13+181 QGLNIHd AJINO παν 1۹ں‎ 
NV Td Z-A NI 3008 Z23 NI 
J'IONV SSdulS NIW SSJUILS XWW ZA ሪሬ እእ ΧΧ ull 
S L&L ٤٢۹ 3 دل‎ 0 7d HW O 9 LNIOd WAG WN 
SNIVHLSOICLOL S459 43015 STGvVL 55.3815 ګل کالکالل‎ 


` 


WALSAS ۱1ا 1۷8010 108000 ا ۷ ل3‎ QGHunsVdW duV SdSSduls 


) 5545 SNV'ld) 
(ΝΩΩΝΙ1ΝΟΟ A/Z) l 31 9 ٣ 9 ھھھء‎ ۹ 0020۹ 11] GUN ΝΤ Ὁ Π.Ο 0ٌ 0۹۵٤ 


3 


00*06- 


3 


CTL ISOS O CO IES ο EU IESE ο- CU STE, ος 
004300000 50+000 0 ¿I [ከ 0> ۹۹۹0ی‎ >6 0 00+0000 6 


a2 


(SALVOIGNI OS (SdI) ٧٧ ٥د‎ LNGW3'I4 N3HM (iLNINd AJNO SHY SNIVUIS) 
3NVid Z-AÀ NI  3NVid Z-À NI 


0 4 5538175 XVH ZA Z4 AA XX " HELI 
es 1 ο መስ 1  پ‎ 5 LNIOd 0۱ء‎ 
SNIVUHLS 'IVLOL / Sd3SSquLS SIMUL SSIS EE gE 


WALSAS GLVYNIGYOOD TWHO1D NI CIUNSVAN 3۷7 6765 


(553815 311713 ( 
(ΝΩΩΝΙΙΝΟΟ A/Z) C 4 [1 O 3 5 L N ae د‎ do ΟΙ ሮሮ ΙΤ» ئو‎ 5 ae 6 


( 93d 06'06 ) SLV'Td 21403105100 30 7 NI ۸۹٨ ΖΊΝ/Ο `S VNIGV NOISHSA - ۷٢۷1٢ 463 


2-D PROBLEM 


(Plane Stress 


[905 / 7521ሬ) 


*Also includes output of MIC-MAC program 
showing the laminate stiffness Eio (page 2. 
column w, rows 22-28) 
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1 [- .ری تد بس‎ 8 τ. ὍΝ Piet TAER 
2 |_[tħeta) |; -82.5; 825; 825 “82.5! [repeat]: ከ | ከይፊ-3 EES oek] 
41ቾ/20ህ9] ; [119/10 το S. 9. 91-56 =- οἱ 0.00 
5 | 8/8511 —_ 608. "508 595 "508. ከ፲. - 609; E 1.90 
ER R/degraded 627 0 7 627) RAPE 627 RR ہے‎ 
: : ; R/lim i 4.18: R/lim*t i 4.18 
EN NJ Em oc ki .4sigrna? : «sig?lim : rud ‹519›ህ11 ፡ ሼዋግነከ _ P OEwEL 
لا‎ που το ο ο DD ENDE 0.426: 
10] 2 O25) 250; 10451 1045: 1567: — 1724] 08997 ም 
16 ም 7 707 0: -. o 0! 99 ΠΕΙ በከ 
62]. x NN NENNEN سل ماس‎ (٢ 
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[6 6 τ ΠΠ... 000; 669; 2000: 009: ل‎ M 
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ΠΏ ΝΣ Top icQmost | volf i Em i Efx i Xm | Xfx iEm/Em° 
[19 | Baseline 220! 0.005: 0.70: 340; 259 400: 2142: 0.30 
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9 7 . 
. + . 
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. 9 . 


25 | strength £ AREF! AREF! 5 7 ال تی‎ NN 7 ፡ i መ 


27 ኮህ 6:13 liked to Ply tafe στ _ a 
۲13/0521811 191; 103i ο 717 ቄብ 122: 160] 
[29 | 15001 1500: E 246 K 12 2. 0005 18 
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1 |ኮበለር[[ለ LAMINATE MODUL E MODULE - elastic and gems constants, 


[3 | [Ange] ` CHAUD 9597. CESSAT NE. E u NE 

[4 | [thea] T 825; 825: 825 -825 8፻9ጐኦኑ እ 3| ὋὋ 
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.Bottz* 5 DS 050: 025 000: | EE NE 4 
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[21-12 | — 554 — 554| 554 554. 001 4E03 — PRO 
E: 66 |. 981 931 981 981 . 001  7E05 9. 
15 E. 124 12 124 -ዝሪ4ብ — 000:  3t04 ኑሁ፦ š 
20 | 62506 1 m í 2095. 2195 - ፡ 7 τος σου ٦ 
Em Ἢ τν; o. 00ھ‎ 00:70 CS 
Compl i| 218/68: n ντ ا‎ ሰር... eee 
ER  87[>ه‎ ο. ο. MC DER Ta 
8 DM M NEM ο... e ፡፡፡ 
i io ο... i — ው ሚመ. ۰۴ 0-1 
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| 1 [MIC-MAC/ZIN-PLANE: وہ‎ a - ቦክ im 11 تو‎ 
- [ሠ:] | -82.5፡ 8255: - 825. -82. 5: [repeat]; h,” ; ከደ-3 8959] 
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[464/88ክ ; — 1.000; | 1000. “1,000: 1.000; 1,000: 1000; 1000፤ 1.000 
[22 | Mod-Bsin: 00 0.000; Hot/wet} — 340: 259 40.0: 2145 

| /Q- -iso i atlimit ; û ΠΝ ን : 
stiffness 3 πε es 1 1 O NOME .. 
zo ndr E ᾶ x 3 : 


e 00 4 co e 








. ፡-- - ሙ o 














* - መጨው መዱ rera d s aa 


š 
1 
ات مم سس ےب سس ےت مج‎ S کے‎ ሙና س‎ À — — መመመ መመመ መመመ ። سح وس اا‎ መ 
* 


Phy data linked to Ply Data File: 


1፡98 K G NIU ANE 2. 122 160: 
29 | . 1500! 1500: 40: 246: 68: ጅን 22: 00005 





: : : 
. ጊፆር 
^ 2 P 





‹ጁ-ጅ =ው ጠው > ቂጠኡ= —T —— —— —— 


20| απ . 125: 07 16. በጋር 0.5: 02: 0.9: 


— CUm --- o ፡ —- ممه‎ — s وم‎ 


0.02: 225: 0: 0.6: 0316: 0.004: 0.004: 2000! 


ደው ው ው ር... መጨ- iq ርን መመመ sss ባባ مسا نټ ومو‎ ——— — + 
a.s... 0-9. 


A S9.7:_ 2529 NEMO! 1.56: 0: 0! 0: D: 











τ 5 
سو‎ fg ساس‎ 

















29 | 
307 
E 
32 | 
A NEG sa 
ΕΤ 
38 
39 | 
ETE 


1 ۳7 LAMINATE ποσο MODULE - elastic and hugrothermal constants 


مجه ههه و حه صم مم صي چو پوه EQ uu Ie ο‏ جموه ه و موو مع ومو سه سم وچجمچو ووا re ee መቸር sspuns‏ ده شر و ری شش و c Mahe SEE Cow OSS ota EL‏ 


Ns κας سا ہہ‎ Re و‎ NN 
[4 | [et] :— 900 «900: 750; — 750; [REPT]* ١ IM. 
E ap: O ` T MENTO EUR να ο ΥΓ 


.2Xyad : 30400: — SEO — 0 
7 | 4/34 ز‎ 6 E. ue 56:00 55:00 4 





9 | 1፡5 — imb σσ B وه‎ 935 a د‎ 
[10 | Bottz* | — 075 .- 0505 025 000 i  . : 
EB een "κ. ን 0 DO) ٹ٤‎ : 

ο ο ٹٰٰ‎ 1 , 1, 7570 


موم وو የወወዑው69 ወፍ ዑፁውዑ ውው "ውውው ቸው ዕው “ው ጣዊ ወውው=ወ‏ وه هسوسو ዋወ ፍው 6 999/96 ቅፍቁቹ ወፍ ባውወውውውፁው፡ቹ 8 ወብመመጽድውብው ጣፍ ፍፍ ፍፁዓ‏ وو وو ده هوو وه د ده د د هه ووجو د و د ه وسوو وب ه چو و هو ه سه و و ውውቐዐ ሁፍ 9 ውመወ ቁ-ቁፍፍ-ቁዌ 9 6 ፍፍ 6 ውሞ 94-^09 00909 Pb» 00 9-900 49 ወው‏ 
. 


1311 01 ከኤ - ANA. ON MEER DENEN. 
14] 56ከበ ፤ [0]/] : [0/2 i lQ/z i lola i [ለ] ij [Αα] 
mi 
16 


— 11.1. 1035! — 1035: 1198; 1198 — 00! HO | 
- 22 ! 18181; 18181 16047 16047 0.7]  1E0!  — 
πο. τα P0290 CGS 12 09 δου τρ C 
MT RENE: .ያ ቆ፡‹፡ ευ 0 
19| 61716 i 0.00: — 000: — 436i — 436: . O00! 2105; 9. 
|20| 62-26 | — 1! 0.00; 000 , 5650: 3850 , O02; 2t02 à. — 
۴ |! ፒፒ 05 9) ο. ሸረ!) ቢተጋ 5 ፡..:2.: 
22 | Comp ! [3]ጠ/0ከ : HUNC “9ۃ‎ 7۳ 

ET s cU ፳-02ፒ i 10635. —— | — 
E. مو‎ ες... :. π65..... yT 
EN: - 0 MICAT 0ء‎ . 0۳ 
me. i oco MEE DNE ER 51"... ;....- 
እ ٦1ہ‎ 16 ! 22 1 | &— 2111 ከለዉ کر‎ ۰" - 
28 | 62-26 i -1094;  :!. ο πλ ت0 ےت‎ EN _E,GPa 
Nl Ree Os ee EN E د لېه ر:‎ 
Nonrnechanical strese(Pa) and strain Eee SN ee sS 
| 51۱٨ 01: ΕΟ EG ee -2t-001; -9E-01: RRs _ 
NEA. ιν Εμ Περ CUM EO ۱۱۱۲ی نر‎ 
Eo 1. pn/T i pne ; sig ‘n/t í sign Ze i aphao i betao Í eps/iso — 
54] 1 i  2t-04|  4E*00: 2-04: ΕΕ100᾽ 21417: ΟΤΙ! relstiff _ 
55] 2 1-0-05 -28700 75-05: - 26700 0028! 207 
25| 6 i: signo i epsno i  -2t-05; -6E-01; — -5.713; — -0.152) — specR 
s1| 1 . 8-05  JE04| e/x I -2t-06 00011... Miso 
μαι 2 ! . 2E-03: — ፳-07፲፲[ εὖ i  —8t-04 | -0001 ^ 1  fn(nu) 
2 EN ቢር ው Hod ኛና | Ε΄ 


ο 


ADINA-IN τν 7 ٢ FITE 


THV2D. IN ANALYSIS OF Ντο 7٦ 


++ + 


* 


FILEUNITS LIST=8 LOG=7 ECHO=7 
CONTROL PLOTUNIT=CM 

COLORS BCODE=BLUE 

x 


DATABASE, CREATE 
WORKSTATION DEVICE=0 
ል 


HEAD “THV2D.IN ANALYSIS OF UNISOTROPI TEENIE ο σα. 
Ἂς 


MASTER IDOF=000111 
PRINTOUT VOLUME=MAXIMUM P=STRAINS 


PORIHOLE 011111 ΕΤΕ; 1117470 
پل‎ 


COORDINATES 

ENTRIES ODES Z 
[ 2/08 20, 
ΛΟ 20. 
3 20. a 
ھ4‎ ٤7 ο. 
5 20: no. 
6“ 20 = oe 
7 209 20. 
ΤΘ -20. 
3 ο, ZU. 


x 

ዶ 

PIX. ο 9 

ΤΙ; +0 

A 

MATERIAL. lI ORTHOTIROÓBTIG EA-173.360505098B510.45E05 950 DES I P 
NUAB=0.03 NUAC=0.03 NUBC=0.28 GADZSS ٣ 

* 

57-1۶۰ +9 495 

٣01717102 ٠٦٠٦‏ ت۱ 9/052۰٠ ۱٥٠٦۰۰‏ جا 

٣٣٣٦٢‏ ا 121-5 ا5ط 5-7 © 5 ٢۳ب‏ ط50 تا 

STRESS 13 


7ء 

Ἐν πο ΕΕ Τα BEI TABLE 
1 0 8279 1 
7 03 BS l 


x 

LOADS “ELEMENT 

1. 2 7687 ۱ 
" ; 

6:00۳۰ 2 6 ۶۶ 705 
GSURFACE 3 4 CGC 5S 1> 1٣٣٣٢٥ ت٦۹‎ 
Πε] 4 4 ۶ 

EDATA 

ENTRIES ٣ BEL TABLE 


1 0 02.5 l 
X 


x 
VIEW ID=l XVIEW=1 YVIEW=0 ZVIEW=0 
EGZONE NAME=TWODSOLID 

1⁄2 


98 


FRAME HEADING=UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=TWODSOLID NODES=11 ELEMENT=2 BCODE=ALL VIEW=1 
= 


ADINA 
* 


END 
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COSJPOCU SOT ¿O αν τς O EU ODES O 0-66 0 
ዐዐ "ዐይ፦ م٢١‎ 155605 حل‎ S0*d000£'0 Lcd ο ιο” S0+3000€ ` 0 κ τ 9- 0043000070 6 
Í 
(SsSGSIVOIGNI OS (SdI) 9913 LNHW3TS NSHM (LLNIHSd ۸10 3uV SNIWULLS) 
dNV'Jd Z-À NI dNVId Z-À NI 
STONY SSJALS NIW 553835 ZÀ ZZ AA XX ELI 
SL ND SIN Od +0 O o LNIOd ua W(YN 
SNIVULS JVILOL 7 SASSAULS AIHVL SSHubS INAWS TS 


WILSAS ك3‎ ۷ ۸11081000 IVHOJO NI 019/015 د١٤‎ 38ኛ 53553815 


(SSG2LLS any id) 
(WONNTLNOD A/Z) 5 d TORATI L NANTI a SU OE SON 0 1 ۹۹ ۹٦ S S dedu Ὁ 


( 930 5ε΄06 ) GIVId DIdOULOSINN 0ک‎ NI'üCAHL CIN/0'S ۷۷10٧ ٧01531۵ - 8 4 
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CU 7176171 ٦ εὐ Ὁ £0-389£8^0- £ 0 38/66 0- 
00'06- 0030000 0 50140001 0 ۲٦-3900202۶ ۹۵ 6۶ ج0‎ 00+000 6 
ር 
0ظ"‎ - O2 2/۵ CO IBE O. ευ 48266 0- 
007 06- τ απ τος 0- 60 00 ٦٦۳۰3 01067007 0 حل د )9ل ناځ‎ 0013000070 6 
1 
(SHIVOIGNI OS (541) لد سه لد دل ۷ه‎ د٤‎ NSHM GSLLNIšd AJNO SW9V SNIWULS) 
38 2Z-A NI ANY Td Z A NI 
STONY SSJALS NIW SSJALS XWW ZA 22 AA XX LI 
S ዜ N d. N.O d W O ጋ LNIOd 831 
571۲۷۱۳5 1۷07/53 LS 3187.75 5306 INAWIT TA 


WiLSAS HIVNIG4002 Tv8g80795 11ا‎ 030105۷3۴ dHV SdSSduls 


(SSIS 1٦ 
(ΝΩΩΝΙΙΝΟΟ A/Z) τ ۹41 6 ٣ ETN THM 1. 1.1 ἅ ο 3 8) +۹) ۹۶۵۵ ک‎ 5 3 4 LS 
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3-D PROBLEM 
[904]s 


62 


^ ር PDN AP 7٤0 UN PUT PASL E 
x 


K DIN EIGHT SYMMETRIC LAYER USING 3D. 
x 


5۷۱۳۷۲۷۲۷۷۰۲۶ LIst=8 LOG=7 ECHO=7 
CONTROL PLOTUNIT=CM 

COLORS BCODE=BLUE 

* 


ESTABASE CREATE 
WORKSTATION DEVICE=0 
پل‎ 


HEADING ' THV3D.IN; *^*EIGHT LAYERS SOLUTION 90,90 DEG AAR“ 
* 


EASTER IDOF=000111 

PRINTOUT VOLUME=MAXIMUM IPD=1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED-YES FILE-60 

ἡ 


COORDINATES / ENTRIES NODE X n ረ 
1 2205 220 ο} 
2 DUE 20 ον 
3 E20 -5 D 
4 δν. = 0 
5 -20. που 0 
6 ۳ھ‎ J20 ο 
7 ATE 25 Us 
8 20. 2 9. 
9 20 5 o 
10 20. 0 
Jat pU 20 0.1 
12 20 20 Owl 
1 = S O. 
14 20 2 O. 
i 220 20 ο 
16 20. 20 ο 
ከ7 እ 20 On 
18 ME 20. 0. 
101 D O. ο 
(+0۶ 8 ο 0. 
111) 5: J0. 0 
104 Tos 2.68 Ox 
ΤΟ 66 10 ο 
106 5 10 ο. 
΄ 177 15 10 ο, 
-υ- ee LIE τς 2-272 2ل ل2‎ ንን m E ur = 


ہے መ ፍ=5መመ ο‏ ہے ہمہ “πο πα — — — — — — — — — — — — — — — á å å‏ سن — —  —‏ حا ټی سے سے سیت ہی ہی ہم ہم 


- መመመርመመ ہہ‎  — ω سے حم سے سے سے سے سے ہمہ‎ መመመ መመ s መመመ s s s s s s ሙመ።መመ መመ ሙመመመ። مس‎ 


MATERIAL 1 ORTHOTROPIC 1B1.E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 
GAB=7.17E05 7.17E05 7.17E05 

A 

EGROUP 1 THREEDSOLID MATERIAL=1 TABLES 

GVOLUME 1 2 4 3 5 6 8 7 ELl:1 EL2-1 EL3-4 NODES-27 

GVOLUME 9 10 12 11 13 14 16 15 EL1-1 EL2-1 EL3-4 NODES-27 

AGVOLUME 1 2 4 3 5 6 8 7 ELl-1 EL2-1 EL3=4 NODES=8 

AGVOLUME 9 10 12 11 13 14 16 15 EL1-1 EL2-1 EL3-4 NODES-8 

GSURFACE 11 12 16 15 EL1-1 EL2-4 NO-9 

LINE STRAIGHT 17 18 EL=4 M=] NC=ALL 

AXES-ORTHOTROPIC 1 NODI=101 NODJ=103 NODK=102 

AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 

AXES-ORTHOTROPIC 3 NODI-101 NODJ-105 NODK-102 

AXES-ORTHOTROPIC 4 NODI-101 NODJ2106 NODK-102 


63 


AXES-ORTHOTROPIC 5 NODI-101 NODJ-107 NODK-102 
ASTRESSTABLE 1 12 3 4.5 6 ΕΠ. 

2 8ؤ ۶ؤ 9 

EDATA T ENTRIES ΕΕ ο ΡΕ ο 


O YOO GJ BJ 5‏ 
نبربربربرنٹے بے بے بم 
نر بر بے بے پر پے بے B‏ 


* 


LOADS ELEMENT 
l- Z2 ۹ 
10 
ሰ (2 2 074 


* 


FIRE 123 ΠΡ 
PIXB 25 29/1 ο.“ S 


AFIXB 23 LINES / 1١١ ١۹ ١٢۴ IAS + 10٥ 
A 


* 


EGROUP 2 THREEDSOLID MATERIAL STABLES 
GVOLUME 3. 4 10 9 7—8 14 FELI- ELA LU ٢٢24 ٧٢ ٦ 
AGVOLUME 3 410 9 7 6 14 (ا‎ 1,11 FED2o2 ٤۹ 
AXES-ORTHOTROPIC 1 NODI=101 NODJ=103 NODK=102 
AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 
AXES-ORTHOTROPIC 5 NOD1=101 NODJ=107 NODK=102 
^OLRESSIADBEB i) ) 27990 506 DE 1 25g 
311175353111 NTL 
EDATA / ΕΠΕ و189‎ EL VIE SS 

1 1 


 —‏ بے پر ىپم 


2 1 

3 1 

4 1 

ፉ 

* 

EGZONE NAME-THREEDSOLID / 1 / 2 

FRAME HEADING=UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 


MESH ZONE=THREEDSOLID NODES=01 ELEMENT=1 HIDDEN=DASHED BCODE=ALL 
پو‎ 


ADINA 
A 


END 
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9 0 O 
00+39 0 


9066 0 
00+06 O 


90-21656۲51002 
00736۲0017 0> 


9023902729007 
00 4 ο- 


9 0-0 38 0 
πο ۵ی‎ 


20 O 
To - +4 0 


το ο νο” 
TO ICSE OZ 


90-31700۳0 0s 
10-έ0Υ596”0- 


ZA 


ST 4 ۵ 
0٣-000 


9I-dCESLC 0‏ 
© ص۹ 1ء۳‌0""۳) 


9۲276 0 
060 4 0 


٠ 


Q‏ ق80 
ο 6 0‏ 


9 6 ο 
ο) ሠ ٤0 


9 0 0 
0) O 


90 ሠ 0 
00 0 


ο) O 
06 8 0 


ZX 


(ΝΩΩΝΙΙΝΟ2 α/ε) 


ει οσοι ο 
8) آ00‎ 


۶1۲-660 
60-3190 


70 0 O 
BO 0 


11-3765 ο 
80-d988Sv'O 


blo 3901270 
BO 0 0 


νι αοτοςςρ 
8 0 


Pl Joro Qr 
807-0 OF 


۶:۲3 00150۶ 
80-0717۲080707 


ÀX 


S < N 1 NINO € ο ο Ὁ 


50۰7990625+ 02-> 
ርግ 1ከ1ሥርረ O. 


۲0۰۹٢٣ 
20+ 3ت‎ 59 ۲۶۰ 02 


ευ ΒΕ O 
COTESCLLE ٤ 


6 O 
COTACEULE O. 


Fo. CO 
10 71007 ي7‎ 07 


EO 155» 0Ξ 
το ۶ 0= 


5:073 0159017 :07 
[101 ዓሪ ο- 


፣ እ τος ος 
! ποσο τς O2 


ZZ 


SNIWULS 'IWVLOL / SdSSdulLsS 


I d N0 49 


L N 3 S 1 تک‎ 


CO ασε O 
50316966٥٨ 0 


O‏ :5ه ان 
νυ ۵‏ 


ሬ0--ጋፍረፍ91 O 
9َ 0 


ο τυ 
0 66 0 


20-90 
06 0 


COILS O 
0+ 0 


CO 90 
501356 66 0 


CO 0> 
50+8866 6 ۵ 


(SXIV2IQNI OS (SdI) OW14 LNSWSAIG NGHM ቪ3፤ዘ፲83 ዶባክ0 38፪ SNIVULS) 


AA 


aod 


tO-JPPSCP COS 
2۰ +٦ τς 


٢7 -‏ 0-11۷ع 
COTAVELVE 0-7 τς‏ 


εὐοτα νο νο 
COTAVE LIS ο- Lé 


۲٣0۰۹٢۲۶۲ 0۶٣ 


"0 9 6 0- De 
τυ ορ ος O2 
TOASOE O= τό 
ου ο ο ος ος 
0717101747 τε 


CEO TOVE O 


1 ተር! በ በር 
EO που O. 
103326001 OZ ፲፻ 


2 ALI 
LNIOd 


321۹71 554819 


5 ٢٧١ 0 1 4ل‎ τν τρ 1 v 9 


038001 
ګل کا الال 


WALSAS له‎ ٧٧٢108003٧ 1۷۹۹010 3HL NI GGLV'nO'WO ddV SdSSduls 


S S 3 ۸ L S 
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30-06 0O 
("060٤0 


2 0 
103215656 0 


90-10567107 
۲97915076602 


9070960 0667 0-> 
00706960607 


ZA 


ατα 0 
ء000‎ 8 O 


90-36 ٤0 
τι πο 4 ٤ 


ረብ. SET 0 
Ti ٹیپ‎ Ὁ 


Ll 316 O7 
τι 6 0٦ 


ZX 


(ΝΩΩΝΙΙΝΟ2 0በ7/ር) 


¥¥¥ OAC 06°06 


0 ۰ء “0 
0 80-4 


ህ1>=.12፡1 0 
00-3257 0 


7)6 ۶ 
80 ۵ 


7139 O. 
ከን ጋ ους ο- 


AX 


S. uL ٢ ο ν ο 5 
SNIV.ILS IVLOL / 51575 


č dt 0 dg» 


۱1۵1041105 50۸۷۰1 LHOIS +4» 


ερ αν ο- 
603017 σι ٢ 


۲09 O. 
2 4 ΤΡ 


CO ICS O. 
εὐ ΠΤ 


EU saeco. Os 
601200571 07 


Zá 


LN 3 H S T d 


το O 
“ሀ O 


) 0 
0۵ 


ርህ. 55110 
و‎ ۷ 0 


ευ ει 
2 0 


AA 


ከ ας |. 


ت 0 ۸ 


EO SDI O> 
ΙΟ! Τ᾽ 


0 7)6 O2 
TOTS O. 


0-3099 0 
TOILE 


COZ ATI O. 
)01730 27 ٢٢ 


(SGIVDIGNI OS (SdI) 2٧ د‎ LNSWaAIS NAHM QGHLNIHd ۸180 HV SNIVULS) 


ISS 1 11 00003 0 3Η1, ΝΙ 


SNI OPI L V] 


CIN/O°S XNIGV 


1 


Le 


[ር 


1 ር 


GLI 
LNIOd 
له‎ 1٩۷٩ل‎ 95185 


38۷ 
LNAWA TA 


0410 11017 217 ۹5 


ης τ ο 5 5 پ٦‎ 


٧٧0151۱٩۳ - 1, 46 
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3-D PROBLEM 
[9025 / 75215 
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n 8ھ‎ ۹ Τη ΙΒ.“ EST: 
Ας 


A THV3D.IN EIGHT 5 ۷ ۲۱۲۱۲٢۱۶۰ ]م۸‎ ٣٣٦٣ 
ھ‎ 


FILEUNITS LIST-8 L0G-7 ECHO-7 
CONTROL PLOTUNIT=CM 

COLORS BCODE=BLUE 

* 


DATABASE ۳ 
WORKSTATION DEVICE=0 
ھ2‎ 


HEADING ' THV3D.IN: 1ط ٭٭٭‎ 61٦ ×ط ۸۷ط‎ 5 SOEUTION ος ος οσα 
κ 


۱ً) “۹۵٦۶ 


PRINTOUT VOLUME=MAXIMUM IPD=] IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED=YES FILE=60 
x 


COORDINATES ۶۶ ۶۶۹ ٣٤٣ x D s Z 
] - 20, = 20, 0.1 
e 20. 22:10 እ 
3 EL. = 5 0.7 
4 U. s. ο. 
5 - δ 20 ο. 
6 20. 220 0 
71 zu. 25 Ox 
8 20. 5 OF 
9 01. 5 01 
10 20 5 oSI 
11 20 20 0.1 
2 20 20 0 
20۔ دنا‎ 5 0 
14 20: OF 
15 20م‎ 20 DE 
ከር 20. 20 0 
17 یت‎ 20 0 
he 0. 26, Oz 
101 5. o os 
ιο. Qr 0 ο, 
195 3 τι. (E 
104 7 ت2‎ 8 Qc 
105 J 556 10 0. 
106 15. 10 O. 
» 107 1 eu OF 
ہے سے سے سے سے سے سے سے سے ہے ہے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے سے ے سے سے سے ہے سے ہے سے پل‎ 


— ae پس‎  — ہے س‎ መ سج سے سے سے جح‎ መ سے سی‎ á- ω سے سے‎ s s x s s á- á- á- á Á- á- á- á- á- å- á å å- سو سے سے کت پټ‎ m سے‎ 


— — — — سے سے‎ -- - --ω سے‎ መመ ። መመ መመመ መ-=መመ መመመ መይመመ=መመመ سب ہے‎ a ہپ پس‎ να <. «απο 


MATERIAL 1 ORTHOTROPIC 181.E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 
GAB=7.17E05 7.17E05 7.17E05 

* ` 

EGROUP 1 THREEDSOLID MATERIAL-1 TABLES 

GVOLUME 1 2 4 3 5 6 B8 7 ELl-1 EL2-1 EL3-4 NODES-27 

GVOLUME 9 10 12 11 13 14 16 15 EL1-1 EL2-1 EL3-4 NODES-27 

AGVOLUME 1 2 4 3 5 6 8 7 ELl-1 EL2-1 EL3-4 NODES-8 

AGVOLUME 9 10 12 11 13 14 16 15 EL1=1 EL2=1 EL3-4 NODES-8 

GSURFACE l1 12 16 15 EL1-1 EL2=4 NO=9 

LINE STRAIGHT 17 18 EL=4 M=1 NC=ALL 

AXES-ORTHOTROPIC 1 NODI=101 NODJ=103 NODK=102 

AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 

AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 

AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 
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AXES-ORTHOTROPIC 5 NODI-101 NODJ-107 NODK-102 
ΠΕ ΤΙΣ | ! 2 34 5 6 7 8 21 26 27 
STRESSTABLE 1 21 


CEATA / ENTRIES EL TABLE MAXES 


O J QUO PU NF 
بر پر بر بے بر بر ہر‎ 
= Ú) Ü) Z = Ü) Ü) بر‎ 


Ἂ 
BOARDS ELEMENT 
W 0504 
10 
٩۹ ον ϱΕ04 


* 


meee 123 LINES / 17 18 
LE 23 SURFACES / 11 12 16 5 


ንን ን ን ο το πο 12 16 7 1611 / 11 15 15 16 
r : 


x 


σος 2 THREEDSOLID MATERIAL=1 TABLES 
۰۰۰۷۶۲۶۰۹٠ 4 10. 9 %Y 8 14 13 ELl=1 EL2=] EL3=4 NODES=27 
ο. ΙΝΕ 4 10 9 — 8 14 13 EL1=1 EL2=1 EL3=4 NODES=8 
AXES-ORTHOTROPIC 1 NODI=101 NODJ=103 NODK=102 
AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 
EXES-ORTHOTROPIC 5 NODI-101 NODJz107 NODK-102 
ΙΕ 2335678 21 26 27 
ΠΙΕΙ 21 
5 ENTRIES EL TABLE MAXES 

1 1 1 


2 1 3 
3 1 3 
- 1 1 
Ας 
k 
BEZONE NAME=THREEDSOLID / 1 / 2 


Meee HEADING=UPPER XSr=1.0 YSF=0.5 XFHAX-22.5 VFMAX-17.0 


MESH ZONE=THREEDSOLID NODES=01 ELEMENT=1 HIDDEN=DASHED BCODE=ALL 
x 


ADINA 
x 


END 
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“00 S 
۲٣٦٦٦٦ B 


TTT Go 
:4 34 =ኳኻ 1 


5 ہے‎ Qar 


v7 جو‎ Qa 


n3 ር. 


τι DT «αι. 
ELT > #5 اا‎ IT 

Ec 2ه‎ «ης وت‎ 

ως οὐ [111145 1 GHITIN 


xw cJ oec EMINUS). dq He TEES (11: በጠር ጫጫ 1 
Εν HOW zz ΤηΗΡ z NOI GSH HI با ڑژپإلا۔‎ 


T3 


A DEDE DEN ٣٦ 70 7 πε 
پر‎ 


۸ 11131.11 ٢۲)6۲۱۷۰۰ ۵ ۷٢ ۶ 1 د‎ ٢ 
x 


FILEUNITS LIST=8 LOG=7 ECHO=7 
CONTROL PLOTUNIT=CM 

COLORS BCODE=BLUE 

x 


DATABASEQGGLEATE 
WORKSTATION DEVICE=0 
ኣ 


HEADING ' THV3D.IN; ^*^EIGHT LAYERS SOLUTION 90,775 ΠΕΡ. 
A 
MASTER IDOT ου πα 


PRINTOUT VOLUME=MAXIMUM IPD=1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED=YES FILE=60 
A 


COORDINATES / ENTRIES NODE X Y ln 2 

i -20. -20. 0.1E2 

2 20. -20. 0.1E2 

3 -20. -5. 0. 1፳2 

4 20. -5. 0.1E2 

5 -20. -20. 0. 

6 20 -20. E 

7 -20. -5. 0. 

8 208 -5, 0. 

9 -20. 5. 0.1E2 
10 20. 5. 0, 152 
1 1 -20. P 0.1E2 
2 20. 20. 0.1E2 
13 -20. 5. Ge 
14 ο 5. 0. 

.0 .20 .20- 15 
.0 .20 2 16 
152 ,0 20 .0 7 1 
.0 .20 .0 18 
.0 .0 .5 101 
OF 10 OF‏ 102 
.0 .10 .5 103 
.0 2.68 .15 104 
ο.‏ 10 7.68 105 
.0 10 15 106 
ο,‏ 10- 15 107 
کے نت سه هه که که که سے سے سے سے سے جے سے سے لم سے سے سل کے لل کل لل لل μα‏ سے سے نل شک ہے ہے ہے ہے نہ ہے ہے ہے ہے نے کے کے ہے ہے ہے ہے ہے sm‏ ہے -== = ዶ‏ 


— — — — ልመሪረሪረመ፡ጭፋጨመመ መሪ سے سے سے س س س‎ <ፍፕ መመመ መሪ መመል αμ.  — — س‎ oe መመመመመልሚ መመልርመልርመልመሪመሪመ<መል መ 


ων‏ س س س س س سے سے መመረመረመመመልመ። መል á-‏ بحم س 


MATERIAL 1 ORTHOTROPIC 181.E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 
GAB=7.17E05 7.17E05 7.17E05 

x 

EGROUP 1 THREEDSOLID MATERIAL-1 TABLES 

GVOLUME 1 2 4 3 5 6 8 7 ELl-1 EL2-1 EL3-4 NODES-27 

GVOLUME 9 10 12 11 13 14 16 15 EL1-1 EL2-1 EL3-4 NODES-27 

AGVOLUME 1 2 4 3 5 6 8 7 811-1 EL2-1 EL3-4 NODES-8 

AGVOLUME 9 10 12 11 13 14 16 15 EL1-1 EL2-1 EL3-4 NODES-8 

GSURFACE 11 12 16 15 EL1-1 EL2-4 NO=9 

LINE STRAIGHT 17 18 EL-4 M-1 NC-ALL 

AXES-ORTHOTROPIC 1 NODI-101 NODJ-103 NODK-102 

AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 

AXES-ORTHOTROPIC 3 NODI-101 NODJ-105 NODK-102 

AXES-ORTHOTROPIC 4 NODI-101 NODJ-106 NODK-102 
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EECCORTHOWROPIC-5 NODIsS101 NODJ-107 NODK-102 
E M ESSIABEE 1 1 2 ጋ 4 56 7 821 26 27 
EDBESSIABDE 1 21 

ETE 7 ENTRIES EL TABLE MAXES 


CO -J O^ n + CU NJ ፦” 
I2 — => 
F2 UJ QJ F2 r2 owe 


< 
LOADS ELEMENT 
1 11د‎ οσα 
TO 
a ΠΠ. οφ 


x 

Eu 123 LINES / 17 18 

۶۰۰۰۰۰۶ 3 SURFACES / 11 12 16 15 

1280525 LINES / 12 160 /Mo 1l / 1l 15 7 15 16, 
1 | 


* 


ο 2 THREEDSOLID MATERIAL=1 TABLES 

Somer: 5 410 9 7 B 14 13 EDlsl EL2-1 EL324 NODES-27 
K C lh 3 4 10 9 7 B 14 13 EL1-1 ELZ-1 EL3-4 NODES-B 
ο 0۸۶٢٣٢ NODI=101 NODJ=103 NODK=102 

fees ORLHOTROPIC NODI-101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC NOD1I=101 NODJ=105 NODK=102 

መምህ ٗ 109 NODI=101 NODJ=106 NODK=102 

ሙከ ክ110ሀ18ሀጊር NODI=101 NODJ=107 NODK=102 
ESUUESSIABDE 1 1 7۰902 677 

ΓΤ OOSTABLE 1 21 

Cr ENTRIES EL TABLE MAXES 


NJ) Cmn ጭ ርቂ») ኩጋ ት” 


1 1 1 
2 1 3 
3 1 3 
4 J | 


k 

* 

Μο ΟΝΕ NAMESTHREEDSOLID / 1 / 2 

ο. ΠΕΑΡΙΝΟΞΟΡΕΕΙ XSF-1.0 VSFz0.5 XFMAX-22.5 YFHMAX-17.0 


MESH ZONE=THREEDSOLID NODES=01 ELEMENT=1 HIDDEN-DASHED BCODE-ALL 
* 


ADINA 
* 


END 
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PROGRAM CLAREA 


አደ አል አ አአአ አልል ልል ል ል አእ ዴል አል s پل‎ ἐς ἐς ἐς Fx GA GA پلر پل پر پلر پلر پلر پلر پلر‎ Ἂς پر پر پر پو پل پل پلر پلر پلر‎ 


O OOOO 


* 
* 


50 
۶ 


SIFET 


& 


THIS PROGRAM WHICH READS THE NODAL POINT COORDINATE, AS WELL READS 
THE NODES DEFINING THE ELEMENT AREA, AND READS THE ANGLE OF ELEMENT 
OL TENTATION T FROM FILES CREATED FROM THE MAIN ADINA OUTPUT FILE. 

IN ADDITION THIS PROCRAM CALCULATES THE MINIMUM AREA. 


ο ከ ከ ከ ሺህን THE PROGRAM VERIFY THE # OF ELEMENT(NEL), IHE #0F NODES 
(NNOD). 


INTEGER IEL,NEL,NNOD,INOD 
PARAMETER (NEL=12,NNOD=196) 
INTEGER INODE(NEL,4),NODE(4) 


REAL COOR(NNOD,20) , AREAL,AREA2, AREA ( NEL) 
REAL 11٢ 1 ΡΕ ο 2 
REAL ΝΤ 


GI C TERA10 FNAME],FNAME2 ,FNOUT 


٣۰۰٠٠۰٠٠۶ ο ο ENTER FILES NAMES FROM ADINA OUTPUT FILE?” 
τ: NODAL POINT DATA FN ?" 

meme (>, (A) J ENAMEL 

OPEN (UNIT=16,FILE=FNAME] ,STATUS=’OLD’ ) 

ο ብ፡ το απ) “NODAL POINT DATA ΕΝ : ',FNAMEI 

5 NODAL POINT DATA FN : ',FNAMELI 

ITE CS A) ELEMENTS NODES FN ?' 

TED CS, (A)') FNAME2 

Swen CUNIT=15,FILE=FNAMEZ,SIATUS="OLD ) 

ከ በ] “ELEMENTS NODES ἘΝ : '",FNAME2 
IESE) “ELEMENTS NODES FN :—-',BENAME2 

ο ο ο) “ENTER THE NAME OF ELEMENTS AREAS FILE ?’ 
PEROS,” (A) ) FNOUT 

OPEN (UNIT=13,F1ILE=FNOUT,STATUS= ‘NEW’ ) 

1:1 ۱5 و‎ AREAS OF ELEMENTS FILET: '",FNOUT 

NE KIS) AREAS OF ELEMENTS FILE : ` ,FNOUT 


PNG X,Y COORDINATES FROM NODAL POINT DATA FILE 


ENDS o,30) “NNOD ," XCOOR ,'YCOOR- 
ከከ (3,30) "NNOD , XGOUR ,'YCOOUR 
DOSO I=1,NNOD 
|| 1 ብሥዉር በስረ 013 114 ١١ 16٣ 11١ لا للم‎ τες το دل‎ 
Ip n bi ስኤ ከ ኔን ንቄ) ከ ጅ ተጋ ከ ን. ከከ]. O) NL  IT1D 1115 
٣۰۱٠۰ ٠٢٠٢ 560100 ,ل1‎ 1( 0٩ 2) 
11 το 3,40 ር 11... 1 2) 
CONTINUE 


EE UBHINOCELEMENT'S NODES FROM ELEM.NODES FN. 


እከክ ”። 0). ከ11: "NODEl','NODE2 . NODE3' ,'"NODE4" 

ΜΜ; το NEL ΝΟΕ, Νο , NODES’, NODE4 ' 
AMIN-2000.0 

DO 100 IEL=1,NEL 

مل هل1 )لال لء لد 1م 9 وګ 18 16.17 15 14۱ 1 12 ο ο ο ο. TI‏ 


& er IEL 2) INP (JEL 3) BHODEXIEDS4)STI6,I17,118,119 


ον ο IEDOTNODEOPER J) INODB( IED 2) INODS 07ا15‎ 7 


& INODE(IEL,4) 


O TT ፓ ፓ ን ቁ ፤፡፡ NODE CIEC ከከ ከ ፦ዙዚ1ኮ! ሪ),.111(0)11.(ፅ 121,277 


& INODE(IEL,4) 


HOBECL)-—-— INODECIEL;T) 
ΕΟΟ ΤΕ INODBUDED; 2) 


OS 


NODE(3) = ء٦۲٢٣‎ 
NODE(4) = 110115) ۱١ ٤-٢ 
k 


* DIFINING THE X,Y COORD TTO NODE (IDS ΠΠ. 


WRIIE (5,302 NNOD # XCOORI T ۶٥٣ 
WRITE (3,30) NNOD 001 S 56 E 
DO 80 J=1,4 
DO 70 J=1,NNOD 
IF (CI. EOSTNODE (IERE HEN 
XJ) = ۰۳ 
Yi) = COOR I2) 
ነ።11፡ “4ء‎ J 7+: 
WRITE (3,40 KOJI O) 
GO TO 80 
ΕΠ 
70 CONTINUE 
80 CONTINUE 
x 


^ ELEMEITS Απ σα ο 577 


Bl s η 
B2 Επ .) 
Cl =a = eres} 
CZ = الم‎ 2 Oo) 
ከከ ን 1. ው رر‎ 
88 = ه ه هب‎ 731 
CCI c p] (4) 
CC? = ιο رر‎ 


AREA] = ΑΒ ΚΕ 2۸81 Εν, ; 
AREA2 = 4 ٦ 
AREA(IEL) = AREA] +AREA2 
AEL-AREA(IEL) 

AMIN - MIN(AMIN,AEL) 
WRITE )5,90( TIE 11311 ΠΜ 
WRITE (3,90) 4 DTED ΕΕ ΠΡΟ ΠΕ 

100 CONTINUE 

[|] 110 1=1] NET 

WRITE (13,99) 1,ARERCL) ΠΠ 

WRITE (5,95X9T. ARES) ۱ھ‎ 
30 CONTINUE 


2.0 


8 FORMAT (750A7, 4X) ADI — E 

20 FORMAT VIS aA 

30 FORMAT (/30A7-4X) 44990 559 

40 FORNAL (ZT AA ΕΤ 

90 FORMAL (//1ጸ. አከ. ο ٨ ٢ الا لل‎ ٥٦ 
95 ۶۲۰۰ ۶ πι) ات‎ 


SIE 
END 
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x ch e ehh AG e f A e لو الو قر‎ c k X À X A k K X A 27 


PROGRAM TRANS3 


ἐς ς Pe Fx P PI xxx PII Ix AAA A A AG X‏ لو پو لر پر پو پو چو 


C 
6 
8 
8 
C 
C 
C 
C 
C 
8 
C 
8 
C 
8 
C 
τ 
C 
C 
x 


* 


ከከ ας ος ANGLE (THETA) FOR EA. ELEMENT 


— aw aw aw ae ar = ae a -ωω. - - --- -- -- H  — پټ لس‎ መመመ መመ መ ټس‎ 
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5101115١ THE STRESSES OUTPUT FROM ADINA,‏ 07۸د1110 
۰۱٦۷۰ ን” GLOBAL STRUCTURE COORDINATE SYSTEM, TO‏ ۴۰۳ 2۰۰۴۲۶۰۰۳ 

Piper iNet PAL MATERIAL COORDINATE. AND THEN, CALCULATES THE 
Fi T PILI T PER STRESS COMPONENT AND PER ELEMENT AND FINALLY 
F CIO τες INE TOTAL RELIABILITY. ALL THAT BY CALLING THE 
SUBROUTINE ο... 


NEL ዙ OF TOTAL ELEMENTS. 

οἱ, S2 CIHE SIRES IN GEOBAL COORDINATE SYSTEM. 

THETA ۰۲۷۴0۶۰۰۸۹٣٤٣٦ 

ο ο TS SCG SIRE STRESS IN THE MATERIAL PRINCIPAL COORDINATE. 


۳ + ole REEEABILITY PER STRESS COMPONENT. 
Revie, = RIG; Ine TOTAL RECGIABILTY OF STRUCTURE PER COPMPONENT. 


RE ٠۱۰۰٢٢ ۶ ELEMENT 
RT GOEHESPODADLOSRHEDIABEBEDIYS 
ALF IHE SHAPE FUNCTION PARAMER FOR WEIBULL 


BEFORE RUNNING INE PROGRAM VERIFY 11 THE PROGRAM THE ANGLE THETA 


EUT AND TRE # OF ELEMENT. AND ALSO VERIFY ΤΗΕ MATERIAE PROEPRITES 
GIVEN IN-SUBROUTINE BETA. 


ELLE πο I 7S IS RR ININ IS 7s PR IR PDS PY Dus کو و ا ا‎ PR PAIS E SOS USA A پل‎ PRIN ICTS 


INTEGER NEL 

ΙΕ CNBESIAZ.PI- 314515927 ) 

PARAMETER (ROWS=100,COLS=6,COLS1=7) 

BNTEGER ITHETA(UNEL) 
EEREABONSORONS,COLS),B(RONS,CODS),CCORONS,CODSI) 

DE CE NEL SI (NELTS ANEL ony (NED none (NEL) ,oYZ( NEL) 
Rew, M, (Het ACNEL), ITHETAF (NEL) ΔΕ 

BEUSRACTERASZ20 ENAMEENOUT,ENAMEL,FNAMEZ 


|| ርን ENTERCLINPUT STRESS ΕΝ 2” 
ΓΕ», (A)') FNAHE 

ο πμ ο πο) -INPUT-STRESS FN = ',FNAME 

Ore UNTT=15,FILE=FNAME,STATUS= "OLD ) 

፣ ክክ!" ን ሚሜ ተክ! INPUT ELEM. ANGLE FN?” 
PESOS, TAI TY FNAME? 

 -ፆ ከ ከ ከ ከህ. INPUIL ELEN. ANGLE TN = ;ENAME2 
SPENCOOUNIT-16,FILESEFNAME2,STATUS- OLD") 
DONDE, ^) ' ENTER ΠΕ OF ELEMENTS TS 


ECS, CA) X FPNAMEI 

۰۴۰٠۳٣٠٠11417207 PIBESENAME1, SIATUS="ODLD' ) 

DERE (5,4)  ENPER OUTPUT FN کم‎ 

PELLIS: (A) ) PNOUT 

5٦۰۱۱٠۱۱01۷10213 ,۲۲1866 ۶۲00۰ STATUS- ' NER' ) 

VESEES (13,^) ' INPUT 2۹۹ ENAME 

13.5) “INPUT ELEM. ANGLE FN = ',FNAMEZ 
|| τν ls.) “INPUT PAREAS OF ELEMENTS FN = ENAMEI 
pesce i n ο ουτε ο 7 


9 
ያ 


D0 20- I-1,NEL 
E SNO د‎ 11 ኬኩ ን TI ΕΚΤ IB ASAIO, 


& ΙΙ; و 11 ط‎ TL ITO 153 
መ መከ ITT ፣ ፣!). 1.1) THETACT) = 9020 
νι ΠΗΡΕ I) E0. 3) PHETA I T)us 0ھ‎ 


ΠΕ; 25ο) Ν,ΤΠΕΊΛΙΡ 
TRETAT) = THELIAL PIT180;.0 
CONTINUE 
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yk 


WRITE (1 558) 00175 الاو‎ ٠ 
Ἂ 


ATRANSFORMATION OF STRESSES τι 111) 17۹ CSREES 


WRITE (13,11) 'NEL',' SA , SY “ΤΙ STA ;; 
DOSS T- INEL 

READ (15,4) ΤΕΝΟΡΟ ΡΕ) 3) ١٢ لل 52 (1 1 2د‎ 7 ee 
THETA 1) 5 1٣٣۵ 17-168 018 

WRITE (13,61) ΤΕΕ Νο τὸ 941  + 1 1د‎ ου... 

& THELIAL) 


L:COS(THETAF(I)) 
M=SIN(THETAF (1) ) 
A(I,1)- SAUI)A^CD552) +۲٣ οι ١۹ Oy Cis 
AC(l,2)= SACL) 7۸۸)٭‎ 2(٣ ص۰)‎ ٦هی‎ 1 ١ 
ΞΘ 
A(I,6)= -SAC1)*EAM + Sy( 1) psl a R ٣ 
و175)ھ‎ (>7 ΙΕ 
Aly ቀሪ! (١ 
25 CONTINUE 
WRITE (13,190) የ011, 51 52 353 1337 ترا‎ 
DO 80 ۳ی‎ 
WRITE (13,60) 1,A(1,1),AC1,2) δι 5) ΔΡ O n MED 
BO CONTINUE 
Ας 


ASUBROUTINES FOR CALCULATING THE RELIABILITY 


CALL ΒΕ د‎ BEL) 

WRITE (57%) 1117 ۸ د‎ 
REDON OSS EEE 

WRITE CATS RITE ን 
CALE HEPA CAGE NEC ALE) 
CALL RECIMA B RED) 

CHEOL ECRE(B NEE) 


10 FORMAT (77,110 ο αν O D 
TA FORMAT ZI, ISAL ۸ EU ٠ لا‎ 
15 FORMAT TIR, TET 6. 

25 FORMAT (IA, IAEA TI ο. σι 

60 FOR UR 1707; 35 ۶ءء‎ 7۶۵+ 

61 FORMAT (7,45 1676EII19 0 p ال‎ 


) 
) 


51 ٢ 

END 
ዶተ+4+4++ተተተተተተተ+ተተተተተተተተተተተተተተተተተተተ+ተተተ+ተተተተተተተተተተተተተተተተተ+ተተተተተተተ+ 
A hok cho هو او پو پو لو پو‎ A A hh لو‎ 


SUBROUTINE BETA(A,B,NEL) 

پلر پلر پمر Ἂς Ἂς ας κ Ἂς‏ لو لو لو لو ας ας Ἂς Ἂς Ἂς Ἂς ας‏ لو چو Ἂς‏ لو k Ἂς ας‏ لر چو چاو لو 
PARAMETER (ROWS=100,COLS=6)‏ 
REAL*8 A(ROWS,COLS) ,B(ROWS,COLS)‏ 
REAL (GOLDS) KIC ٤‏ 

< 


* PLY STRESS5 - STRENGEHTRATIONSADBCUDLATIS 
WRITE (13,10) 'ከሂጄሂፒ" SR1', 59۹۹+ ++ +111" 
DO 90 I=1,NEL 

C STRENGTH PROPERTIES FOR CRAPHITE ΕΘ IE 


ቀቁ ፀፃ ου 9 می‎ ቀ ው» ቀ 9 9» ፍ» ፍው ፍው ቁ ቁ ቁ ቁ ቁ ፀፅ ቁ ቁ ቁ ቀ» ቀ ፍ ው ው ቁ ፍፍ ቁ ቁ ፍ» ው» ፍው ቁ « © ው» ቁ 9 ቁ 94 4 ®» 


80 
90 
60 


ΓΙ; 554 

X(2) = 0.4E4 

E SM 

X(6) = 0.68E4 

X(4) = X(6) 

SJ = CO) 

ICE = 152054 

X2C = -2.46E4 

ο ο አራ. 

Peel, Lest 0) KOL >=. 
ΠΕ, ΤΙ 2.0) X (2) = τος 
Πο ο ο ο) XS) ር ጋር 


DO 60 J=1,COLS 
B(T J) = ABS(ACI,J) πα) 
CONTINUE 
ome heels. 60) 1700 1,.))7D00(1 2) B 13) 81:6, 1,5,9) .( 
CONTINUE 


ΓΤ ርን CIR EIU OAT, OR) ΙΧ SOC) 
Εμ Αα τσ ο 1.) 

RETURN 

END 
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APPENDIX P 
APPLICATIONS 


We are again considering the simply supported anisotropic square 
plate with uniform tensile load, described in the former appendices. 

Here, the meshing used is modeling a conventional repair of a 
damaged structure. We are thereby calculating the reliability for differ- 
ent configurations. 

Using the THP3D4 program, the verification problem (THV3D) is 
solved with different mesh sizes (i.e., different element sizes). Com- 
paring the results for the total reliability, we may conclude that the 
ΠΠ are in good agreement with each other (RT(THV3D) s 0.807, 
RrtHP3D4) = 0.824). 

The THM3D4 program analyzes the parent structure with a hole. 
uri esulting total reliability is O (Ry = 0.0). In other words. the 
structure will probably fail due to the given applied load. 

The THM3D4 program simulates a repair configuration using a 
weaker material (E-glass/epoxy instead of graphite/epoxy used for the 
original parent structure). The reliability decreases in comparison 
with the original undamaged structure, Ry = 0.1. 

To improve the total reliability, in the THA3D4 program we sim- 
ulate a repair configuration by using a different lay-up orientation 
Wl for the patch and, as expected, the repair causes the total 


structure to be more reliable, Rr = 0.87. 


8 


In summary, from all the individual results for the computed reli- 
abilities listed in the programs below, we can interpret the critical 
element (see RE in the results) or the critical stress component (see 
RT1, RT2, etc.), which may represent the fiber or the matrix of the 


composite material. 
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20 IIC 
CYLINDRIC 
CYLINDRIC 
ء٣۴‎ 
ο ΕΡΕ 
CYLINDRIC 
CYDINDRIC 
CYCINDRIC 
CILINDRIE 
CYLINDRIC 
CYLINDRIC 


GVOLUME 7 4 18 17 
GUOBUME 4 1 11010 
ο ο. T 2 1211 
GYVCEUME 2313 12 
601/1۳ 3 6 14 13 
GVOLUME 6 9 15 14 
GVOLUME 9 8 16 15 
0001115 7756 


COUP PAGE +89 + ۲ی‎ 
7ئ‎ 4 6+ 64 8 


28 
30 
31 
97 
33 
34 
35 
36 
37 
38 


EL=1 
EL=1 
EL=1 
EL-1 
EL-1 
EL=1 
EL=1 
EL = 1 
EL-1 
EL=1 
EL=1 
ELI 
EL-1 
EL=1 
EL-1 
EL=1 


MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 


.270 
ےرڈ 
.360 


NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 


MATERIAL=1 TABLES 


24 


EL1-21 
EL1-1 
1 
EL1=1 
EL1=1 
EL1=1 
ELl=1 
EL1=1 


EL2=1 
ECZ 
EL2 1 
EL2=1 
EL2=1 
EL221 
EL2=1 
EL2=1 


OC aS‏ 0 0 © هته هت یي 


EL2=4 NO=9 
EL2=4 NO=9 


AXES ۶۵1۳۷ 
AXES-ORTHOTROSPIC 2 
AXES-ORIHOTRORIC 3 
AXES-ORIBHOTIROPIC.4 


NODJ=103 
NODJ=104 
NODJ=105 
NODJ=106 


NODK= 
NODK= 
NODK= 
NODK= 


AXES ος αλα 5 
ASTRESSBABEE 1 1 2 
SIRES SLAB Ee ae 


۲717 EL 
1 
81817 TO 
29 
- 2 
STEP 4 TO 
30 
2 
SLIEP 4 TO 
31 
۱ 4 
SLEP 4 Το 
97 


NODI=101 
3۰4.6 TB 7S 


NODJ=107 
2 


TABLE MAXES 


1 


l 
1 


سم سم 


= = 


1 


(ል) GJ UJ Fe 


Mm UJ 


102 


NODK = 


— 


EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 


102 
102 
TOZ 
102 
102 


ΝΟΡΕΡΞΟ 2 
NODES=27 
NODES=27 
NODES=27 
NODES= 27 
NODES=27 
NODES= 27 
NODES=27 


LOADS ELEMENT 


* 


TO 


9 


16 


2 


πω ΟΕ 4 


22 94 


EIXB I2: HINES / B 28 
EIS SURTACES ۷٠٦۷٣6 6 7 8 
پو‎ 


* 


EGROUP 2 THREEDSOLID MATERIAL-s1 TABLES 
ΜΠΕ 11 12 55 31 32 25 25 ELJ=1 EL2Z=1 EL3=4 NODES=27 
ΜΕ 12135 5 ጋጋ 33 25 25 ELl=1 EL2=i EL3=4 NODES=27 
BvOGUME 13 14 5.5 33 34 25 25 EL1-1 EL2-1 EL3-4 NODES-27 
GVOLUME 14 15 5 5 34 35 25 25 EL1=1 EL2=1 EL3=4 NODES=27 
۰.۰۱٣15 16 5 5 35 36 25 25 EL1=1 EL2=1 EL3=4 NODES=27 
۳۰۰۰۷ 16 17 5 5 36 37 25 25 ELIsS1 ELZ=1 EL3=4 NODES=27 
GVOLUME 17 18 5 5 37 38 25 25 EL1l=1 EL2=1 EL3=4 NODES=27 
۳۰۰۰٠۱٦ 10 11 5 5 30 31 25 25 ELI=1 ELZ=1 EL3=4 NODES=27 
AXES-ORTHOTROPIC 1 NODI=101 NODJ=103 NODK=102 
AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 

NODI=101 NODJ=107 NODK=102 


Mees -ORTHOTROPIC 5 
2 


۲۰۲۴۲۶۶5۲۸8 ت‎ 1 1 ΤΟ ΕΙ 


51 ከ ከኮ ንተ! ከ11 ርእ ር 1 21 
ο ο / ENTRIES EL TABLE MAXES 
1 1 1 
SIEP 4 TO 
29 1 1 
2 1 3 
ΞΕ 4 TO 
30 l 3 
3 l D 
STEP TO 
31 1 3 
a 1 1 
BIEP 4 TO 
32 l l 


Ἂ 

A 

ERSONE NAMESPARENT / 1 / 2 

Meni HEADING-UPPER 4SF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=PARENT NODES=0 ELEMENT=1 HIDDEN=DASHED BCODE=ALL 
Zarek HEADING=UPPER XSF=1.0 YSF=0.5 KFMAX=22.5 YFMAK=17.0 
*MESH ZONE-PATCH NODES-0 ELEMENT-1 HIDDEN-DASHED BCODE-ALL 
& 


ADINA 


۸ 
END 


Jo 


Ας Αι Αἱ Pn n eR P P P P e پل پل لر پو لر پر پل پو لو پل پل لو لر لر پل پو پل پلر پلر پلر پر‎ Αν پل پل پل پر پل‎ II AI A A 


PROGRAM TRANS3 


Ας‏ پل پل III PA‏ پل پر پل پل لو پل Ακ ἐς ἐς ἐς‏ پلر لو پو پو پر پر پر پل P xx Ps P‏ لر لو ἐς‏ پل Ἃς‏ پلر پر پو و و لو لو هو 


τ 
C 
p 
8 
C 
C 
2 
C 
Ü 
C 
C 
9 
e 
C 
C 
8 
ë 
5 
x 


አ 


THIS PROGRAM TRANFORMS THE STRESSES OUTPUT FROM ADINA, 
EXPRESSED IN IHE GLOBAL STRUCTURE COORDIN Ir ንዝ! τὸ 

THE PRINCIPAL MATERIAL COORDINATE. AND THEN; CALCULATES THE 
RELIABILTY PER STRESS COMPONENT AND PER ELEMENT AND FINALLY 
EVALUATES THE TOTAL RELIABILITY. ALL ΤΗΣ ٢٢٢٢٢ ل٢‎ 
SUBROUTINE PEREI 


NEL # OF TOTAL ELEMENTS. 

SX,SY,... ,SY2 : THE STRESS IN GpOBAL E OO DI I ፡ 

THETA : THE LAYUP ANGLE PER ELEMENT. 

51,82,... S6 : THE STRESS IN THE MATERIAL PRINCIPAL σου. ss 
R1,R2,...,R6 f WHE ×۴٢ ظط‎ ٣٢٢٢١٦+49 ο ٦59 
RT1,RT2,...RT6: THE TOTAL RELIABILITY OF STRUCTURE PER COP P i as 
RE : THE RELIABILITY EER ELEMENT 

RT : THE 0 

ALF THE SHAPE FUNCTION PARAMER FOR WEIBULL 


BEFORE RUNNING ΤΗΕ PROGRAM VERIFY IN THE PROGRAM THE ANGLDE ο... 
INPUT AND THE # OF ELEMENT.AND ALSO VERIFY THE MATERIAL PROEBPBISEES 
GIVEN IN SUBROUTINE BETA. 


پو پل Pe P B IA AAA‏ پر لر لو لو Ἄς‏ و چاو و پو پو لو لو پل پو پو لو پو و لو لو و پل لو ἐς Ας‏ پر پو e x Pe e e e‏ لو لو پل Ἂς‏ لو لو 


INTEGER NEL 

PARAMETER (NEL=64,PI ) 

PARAMETER (ROWS=100,COLS=6,COLS1=7) 

INTEGER ITHETA(NEL) 

REAL*S A( ROWS ,COLS) BC ROWS, ۹  ከሱፕ )ል)ያቻዯፕ ος) 

REAL SA(NEL) ,SYCNEL) /SZUNEDL) SS (NEL) , SXZ2QNER ο NEE 
REAL LM, THETA NEL) 1۱١ ۱۹١ IEE ٢٢ 

CHARACTER*20 FNAME PNOUT, PN AME 5۵۷٥۰ 


WRITE. (5,^*). ENTER INDUISSISESSEDI EM 

READ ሐ ٧٢ ٢۳ 

WRIIE (5,4) “۲۲۱ 6 ۶٦/١) ٰیٰٰ٣‎ 

OPEN (UNIT-15,FILE-FNAME,STATUS-'OLD') 

WRITE (5,^) “ENTER ከር ٣٦٦+7 
READ NS. ο ΕΕ. 

WRITE (5,^) ٠ ۲۷۲۱۷۰۶۷" ٢ F NATE 
OPEN (UNIT=16,FILE=PNAMEZ, SEATUS- OLD ) 

WRITE (5,^) “ENTER AREAS ος Εμ ΡΕ. PIR 
READS 1 

OPEN (UNIT=17,F LLE=— NAMEI +ى,‎ 1 0 

WRITE (5,) ENTER 0۱1 1 7٦7 

READ CS, ህህ ποπ 

OPEN (UNIT-13,FILE-FNOUT,STATUS-'NEW') 

WRITE 15 INPUT ΕΊΤΕ ΘΠ 
WRITE (13,252) INPUT ELEM ANGLE TOS EEL 
WRITE (13,^) ‘INPUT AREAS ۳ ١ ۵۶۳ 
WRITE (13,4) m OUTPUTS ENT 7٣٤ 


AINPUT OF ANGLE (THETA) FOR EAM ELEMENTS. 


20 


00 0 
READ {16,4} N,12,13;,14 16 ΤΠΕ. πο 1:٤ +7 


& m 1233 1 ھ۸‎ ΠΡ 
IF ΠΤΙ .Ep. 1) ۲۱۰۶/۱۸) (۸۶٦ 
TP 11 ር 1 m ).( 111 1( 2) 


WRITER S, 20I N THETACT) 
IHETAF (I) = ΤΗΡΕΙ; ΕΙ 
CONTINUE 


104 


* 


۰۰٠3,00 *'.FNOUT 
* 


ΝΕ STRESSES TO PRINCIPAL AXES STRESS 


ΙΙ, ΡΠ. Nel له ته هه‎ wo ome , SY ; ANGLE 
ከ ን I=1,NEL 

D EN ل‎ 115071 5175111 5456 1 SXY(1),52X2(1),SY2(I) 
۰۰۱۰۰٠۷۸۰۰٠٠٠۹ ٭‎ ٢۲۸7٦٤ 

ο τα ο ο ος ο τὸ Ρα ο 1 SAYC οκ τὸ οτι 1), 

& THETA( 1) 


E-COSDDPHEIAFECI) 
DESOINCPHEIAFCIO) 
1 (11.1 ۰م(‎ τι σσ) ZASXY(I) ALAM 
د‎ ο τ. ο ο ο ος ον ο )ودل اح‎ ١٢ 
ΞΔ] 
Pie oO ..፡.)፡) 111...) 1 ተ. CIDA LAA 2 115) 
ΠΗ + SXZ( LIAL 
2۰٠۸۰ کو ال ج2‎ = SALI) A 
25 CONTINUE 
ο Εν ΠΤ U51',*92 4053 So MESS. 54 
ο 6O 151 NEL 
Got eel oO; 1٠1,1۸ د1642 ۸ل 1 ۸ 2 و1‎ 977 ACI 4) 
80 CONTINUE 
< 


^SUBROUTINES FOR CALCULATING THE RELIABILITY ; 


s BEIACA,B,NED) 

ΓΞ 5.1] "IHPUT ABFEA' 
p (5,*) ALF 
Νε; ALA: CAT 
Seu ALFACA, B;NEL, ALF) 
SOL RELI(A,B NEL) 

O CCRELCRE(B REC) 


10 mom (//;143,AI0,6(A7,4X 


MUL BO - 
| moet (//,1X,AL0,7CA7,4%) 


DP 
ΑΕΗ 9 
۳ 
6 la) 


15 r ο το, THETA τὸ) 
25 ο ο ο GUAR, IHEIAC το) 
60 ሽሽ... αν 16, 6EII1.2) 

61 46 ΙΕ 6Ε11.3,54.56.1} 


STOP 

END 
K++++LILLLILTLLTiLLiLiLiLiLiLi+iLiiLii+iiiLiTi++i+iiiLIii+ii++i+++++++++++++++i++ 
V PG FU Ph Fh eoe e e e e ee e e e I پو بر هلو‎ 

EUBBOUTINESBELACA.B.NEL) 
پو پو پو قو پو پو لو امو پو‎ Pe RA e e e Ἂς Fe e Ας Ας 

T PARAMETER (ROWS=100,COLS=6) 
FELE BCRONS COLS), BROWS COLS) 
ΓΕ COLS) AIC; X20 AC 


gees 1د‎ 55 - STRENGTH, RATIO CALCULATION 


TIE CONGO ED ሕቢኪይሱጹሱጅች›ጅ STG SET SRO DRST ۰ 
DO 90 I-1,NEL 
"ΓΤ PROPERTIES FOR GRAPHITE EPOXY, 


ቀ9 ቀቁ ቀቁ υ ቀ ቀ م‎ 9 9 0 9 9 9 oo > o o b» 9 » ቆቅ ቀ ቀ 9 ቀ 9 9 9 9 9 ቆ 9 9 9 9? υ ቆቅ ቆቅ وچ‎ ቆ % ο ο ο 


80 


30 
10 
60 


X(1) = 34 

X(2) تد‎ 0.4 

9 ሪች) 

X(6) = 0,68ጄ45 

AO =>) 

72٢9(3 “πω ιο 

ΤΕ ΞΕ 

X2C = -2.46E4 

2 + τε 

IF 1}; ΟΕ ወባ 
IF ACIZ LIO TOL) የመ... 
IF IACI. LIO 9 9. .- 


DO 80 J=1,COLS 
BCI, J) = ABS(ACL I 
CONTINUE 
WRITE (13,60) L,B(1,1),B01,2) B (1,35) 5 1 n nn s M 
CONTIRUE 
FORMAT (7/2500 oU ATEM b) 
FORMAT ΑΘ 6 2) 
RETURN 
END 


Ας k k k k k k k k k k k k k k k oo AG AI‏ پلو 


SUBROUTINE VALEAUA, BaNEL, ALE) 


Ἄς‏ پل I‏ پر پل پو بل x Ab Ade x ex I‏ یلو 


80 


90 
10 
60 


PARAMETER (ROWS=100,COLS=6) 
REAL#6 A( ROWS, COLS), BC ROWS, COOLS) 
۰" ی۹‎ ۷ 


STRESS = STRENGIO RATIO RISED BY ALES ۸٢ ION 
WRITE (13,10) ከኮ” V2 DIM c 
DO 90 151 NEL 
ΕΕ πα ο ሠ Το ΠΠ 

AMIN-400. 
WRITE (5,95) 11,1 ۹ 77:۶ 
ARITE 03,95) ILE TEE ١١٢ ٢٢ 
0 0 8 

ΑΓ ο) το ο ٢٢١١ ( 

ኣ(1,0ህ] = 10)(1,ህ)፡ዳዯ ሸህረያ 760 

CONTINUE 
WRITE (13,60) 1,ACI, 1) AUI ΑΕ ο ο ο ο. 
CONTINUE 
PORATA A / IT ATO 68 T T CNN e (C ON) 
FORTATA ο SIC GEIL ΠΠ 
FORMAL (AIA ae "تک‎ 2 ۴۰0۳ 2 
RETURH 
ENO 


k k k k k kk k kk k k k k k k k k پل‎ k k X بل‎ kk k Kk k پل بل‎ xA 


ሥ 


SUBROUTINE RELI(A,B,NEL) 


k k k k k k k k k k k k k k k k k K k k k k k k k k k k K‏ لو 


k 
^ PLY 


PARAMETER (ROWS=100,COLS=6) 
REAL*8 ACROWS,COLS),B(ROWS,COLS) 


SLERSS TENSOR RELIABILITY CALCULATION (WEIBULL) 


WRITE SG) 10) NEL “Pile Re ee πο ΠΕ, πα. 
DO 90 I=1,NEL 
DO 80 J=1,COLS 


KOG 





τ. E= 0) J) ~ 
80 CONTINUE 
RITE (13760) L B(1L,1) B(I(T1,2yY B—-(1,3) B II 16) 7 B(—(T1,5Sy; RB(1,4) 
90 CONTINUE 
10 BerHAr (//,1%,A10,6(A7,47%) 7 1X,80(*-"*)) 
60 FORMAT (/,4X,16,6E11.3) 
RETURN 
END 
Ας Ας Ας Ἄς Ας ας Ας Ἄς Ἂς Ας Ας Ἄς Ἄς P PP پل‎ xxx ας Ἂς Ἂν 
SUBROUTINE ELRE(B,NEL) 
Ας ἧς Ας Ας ας P Ἂς Ας ας Ἂς Ας ας Ἂν ας Ας ς x x xxx مر ار پو‎ 
PARAMETER (ROWS=100,COLS=6,COLS1=7) 
REAL*8 B(ROWS,COLS),C(ROWS,COLS1),RT 
* 
*ELEMENTS RELIABILITY CALCUL., AS WELL AS "FIBER" (R1) AND "MATRIX" 
| 8 7 CALCULE TOTAL STRUCTURE RELIABILITY CALCUL. 
η: IO NEL RI ۹ 1 ከ3 ΕΕ’ 
RT =1.0 
C(NEL+1,1)=1.0 
DO 90 I=1,NEL 
NIC COELSI) s 1.0 
የ 50 J-1,COLS 
ο τε 1) 
e CoL Me 051 s GUI. Jg. 
CNEL MONE TELLIN A C(I J) 
80 CONTINUE 
ο RI > C(I οσους] 
ο πο ئ3‎ 66757۷۰11507 Oro οτι ο. οι ο) ο ιτ ισα τ) 
απ ο εντ .ጅ).፣ Z2) ۱ ٢٢ د‎ ιτ. ο ιτ 661 7( 
90 CONTINUE 
ESSE ο TON whe la RIL RIZ PRT: RIG CRIS 


& "ظ0‎ 
ο ο ο ο سا 5و‎ ٢٢ ۱١ ل٢‎ ο) ο ΝΟΕ لل‎ 3); C (NEL+1, 64); 
& τει. 5) ü PL 167 0 


WRITE (5,70) PT 
WRITE (13,70) RT 
10 FORMAT (//,1Z2,h10,7(A7,4X) / 12,95('-')) 
60 FORMAT (/,4%,16,7E11.3) 
65 FORMAT (/,11Z,7E11.3) 
70 FORMAT (//,1Z,'TOTAL STRUCTURE RELIABILITY - ',E11.3 /1X,40('s')) 
RETURN 
END 


| 


INPUT STRESS FN «ος ο ος 

INPUT ELEM. ANGEBSBNOCOTUESDC- ED 

INPUT AREAS UF ELEMENTS TN ΠΕ ARES 
QUTPUL ٣٣ 


108 


NEL SX SY S2 SXY 5፳2 SYZ ANGLE 
1 211 ος .391E*05 18 7E1002 2235E 05 S90SESQUO و‎ ٦ 20 
2 ۹۷ء‎ 20 05 ... 7 . 2364 .849E+00 .604E+01 το. 
3 - 2795103 600ر 5 ھ.‎ (7 .236E404 .849ጄ+00 .604E+01 75. 
á +5 .391E405 .187E+02 4 OSE 0O 622E mO 90 
5 . 278E+03 ግ ከከ Os ΡΕΠΟ . َ ۳ Δ ۳ 90 
6 +۹ TOSE EUS ΕΠΙ .235E104 . 36 1١١٢ .683E*01 T7595 
7 «ο Ετος » 209E+05 ء7٦‎ +2 04 Αθ ος ء۵۳٣٣‎ do 
8 :2 28E S03 ΕΠΕ >7 ΕΠΕ مدق‎ .675E+01 307 
5 SEE ES ΕΒ SSE Ol ጋ ው ት ነፈ “ጋጋሪ 11 .216E*02 gor 

10 .278E*03 +2 ۶ -746E+00 . 235+4 . ΤΙ ھ2‎ ۷ To 

11 OEE OS 96 05 +6 0 . 254 6ھ 7 چ2‎ ٣ 75. 

12 27817 9 5 GSE D 26 ۹ +٦ 7ھ‎ 90; 

P3 -2I BETO . د91‎ ZI SELLO -239E -04 πὶ . 210E. 907 

14 .۔ ۹ :. 03 2-۔‎ 0 2 7 .559E+01 2 ۶ PSR 

15 :278E 403 ο Εναν . 9 0 να. 5Ετοἑ .559E+01 21 0E X02 Ts 

16 5ء‎ s 6 ۳ A1 35209 2:2 909 ۰ 9 ٦۲ “ሪህ... 307 

17 va BIS SOS LES 05 .250FE101 . 235E οἱ SSE TOU -434E+01 gop 

18 ¿27551 03 ΕΠ 2 190Ε10} 235E 04 -967E+00 .438EËE: 01 του 

19 2247102012. ESSE US . 218E+01 ۵۶ .967E*00 -438E ON 295 

20 5۵۹و‎ ን ۶۹ 0ؤ‎ 1111 235E 04 .992E+00 .434E+01 90 

271 .278+03 ISIE OS SISUEUCOA “×2 +0 09 .566E+01 ο 

22 N27 ESOS ο» .163E*02 .235E+04 -550E+00 4572 3E TNI TOR 

23 መሪ! .5 07ھ‎ .2355t04 . 0 ء7۲۷٢‎ 757 

24 275685+03 2ء۸‎ +55 . OP . 236E404 .608E+00 . [1 30 

25 .318E*03 -390E+05 115E+03 .235E+04 [076 .109E+01 3907 

26 2) ETOS . 209+05 70, 982E E02 068 .607E+01 .414E+00 75 





21 
28 
ሪጋ 
30 
31 
E 


14 
5 
16 
17 
18 
13 
20 
2] 
22 
23 


24 


224 3E403 
. 310 3 
+03 
-241E+03 
στου 
20 03 
2015303 
+۹03 
ھ8‎ 1 Ετος 
Eo Gg lES03 
۶03 
1۶975۳+03 
2 2 +3 
۶۱+03 
03 
ο ΕΤΟΣ 
ο ያጋ ህጋ 
2 ۶3ء‎ 
.289 3ه‎ 
2274 E*03 
. 274E*03 
2202E-03 
2 ሰኮ+0. 
.263E*03 
22635403 
۰903 
. 295E«03 
. 264+03 
. 264+3 


. 225503 


2095205 
. ۶ 
098 05 
ο Ετος 
ΟΠ 
SOE 
8 ος 
0 ی۹‎ ን 
ጋ መ Ετος 
391805 
ος ος 
2209E 09 
. 209+5 
ΞΡ 
,+, ۶ 
ο Στου 
:209ET05 
۸.5 
.391E+05 
:209F4 05 
. 209+5 
. 5 
"3915405 
29 9 ος 
3 01. 
115د‎ 
ΠΕΠ 
. 209F+05 
USES 


. 3901+ 5 


:982E*02 
؟؛؟11|5.‎ 3 
72ء‎ 
.S80F102 
1900807 
70+1 87 
۰594400 
.10 7۳1 
.107 1 
. 90 
.646E+00 
.117E+01 
+۷ 
-. 0۵ 
.320Ε00 
.274۹ 00 
. 274 0 
+9 ٤٣١ 
017 
9:9 ٦ 
+ ۹٦ 
.641E+01 
ግ ሪት (ረ. 
. 361۳7 
.36 171 7 
.422E*02 
. 416+2 
356E ۶ 
. شا د‎ 


1162ء 


. 235+4 
ي7‎ ۷ 
۵ 
. 235+4 
. 235+4 
. 398 
22550708 
.235E104 
7428 
, 94 
- 2355104 
. 235+4 
.235E*04 
- 235E+04 
"295810. 
. 235E+04 
. 43504 
:235E*04 
. 235+4 
-2366+04 
. 236E+04 
«2355104 
ገ 04 
.235E*04 
"¿2555104 
۶ھ“‎ ۷۹ 
.235E 104 
.236Ε04 
"2388105 


. 2355+ 4 


109 
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1ος. 
21 
τον 
315: 


116 


C) 2 ©‏ یی یم یی یی می 


AL 


g = =‏ پرب بر بپر = 


Á 


INE CYLINDRIC 1l 
PUE CYLINDRIC 12 
PU F CYLINDRIC 13 
LINE CYLINDRIC 14 
IU CYLINDRIC 15 
CINE CYLINDRIC 16 
CINE CYLINDRIC 17 
mee CYLINDRIC 18 
fee CYLINDRIC 31 
ΠΠ CYLINDRIC 32 
۰۰۰۰۲0۷51111510 93 
LINE CYLINDRIC 34 
LINE CYLINDRIC 35 
ENE CYLINDRIC 36 
NE CYLINDRIC 37 
۳۰۲۷٢۶" .10ذ ۷111107 ت‎ 8 
ዶዴ 


E ROUE 1 THREEDSOL 


18 10. 360° 


30 ιο. oO 
37 IOS ES. 
32 10. πο 


33 του 155; 
34 ιν MEO 
39 10. ο. 
36 .ل1‎ 20 
237 TO. SSIS. 
38 T0 E359: 


10 EL=1 MIDNODES=1 
1l EL-1 MIDNODES-1 
12 EL=1 MIDNODES=1 
13 EL-1 MIDNODES=1 
14 EL-1 MIDNODES-1 
15 EL=1 MIDNODES=1 
16 8121-؛‎ MIDNODES =1 
17 EL=1 MIDNODES=1 
30 ELI ٣٦ 
31 EL=1 MIDNODES=1 
32 EL=1 MIDNODES=1 
33 EL-1 MIDNODES-1 
34 EL=1 MIDNODES=1 
35 EL=1 MIDNODES=1 
36 EL-1 MIDNODES-1 
37 Eb-1 MIDNODES-1 


NC=A 
NC=A 


ID MATERIAL] TABLES 


27 24 38 37 EL1=l 
"ھ2‎ ۹۶ 


ሽን =>. 
EL2=1 


νο ο ο ته ہے‎ Ooo 


H|- 


EL3=4 
EL3=4 


C OCUME 7 4 18 17 
ΟΠΕ 4 1 11 10 
BSBUME 12 12 11 
ο ኾኑ ሪ 3 13 12 
GVOLUME 3 6 14 13 
GVOLUME 6 9 15 14 
ο μυ 9 8 16 15 
S OLUME 8 7 17 16 


ESURMPACL 7 8 28 ሪያ 
ESNEEACE 8 9 29 28 
ZESSS-ORIHOIROPIC 1 
"۰۹ 01۲01٣٢ 
۶۰۱۶۰۰۰ 1.. 


AXES-ORTHOTROPIC 


2 
3 
ES ORIHOTIROPIC 4 
3 
Z 


A I F ESSTABLE 1 1 


STRESSTABLE 1 21 
EDATA / ENTRIES EL 
1 
STEP 4 TO 
29 
RC 2 
STEP 4 TO 
30 
3 
STEP 4 TO 
31 
4 
STEP 4 TO 
32 


20057 21257: 88 8 9 


22 
23 
26 
29 
28 


ΕΞ 


32 ٢٢ E 
33 11-1 


36 ٤ 
EL2=4 NO= 
EL2=4 ΝΟΞ 


NODI=101 NODJ=103 
NODI=101 NODJ=104 


NODI=101 


NODI=101 NODJ=106 
NODI=101 NODJ=107 
ο ο - 6 


TABLE MAXES 


1 1 
1 1 
1 m 
1 3 
l 3 
I 3 
1 1 
1 1 


NODJ=105 


EL2=1 EL3=4 
EL2=1 EL3=4 
EL2=1 EL3=4 
EL2=1 EL3=4 
EL2-1 EL3-4 
EL2=1 EL3=4 
9 
9 
NODK=102 
NODK 2102 
NODK=102 
NODK=102 
NODK=102 
27 


NODES 727 
NODES -27 
NODES 227 
NODES 7-27 
NODES =27 
NODES=27 
NODES =27 
NODES =27 


LOADS ELEMENT 
Y 2 3290-04 
10 
o eS 


۸ 


FIXB 123 LINES TS BB 28 
FIXB 23 SUREACES:'/ 2708.28 27 BTS 
* 


* 
EGZONE NAME-PARENT / 1 


FRAME HEADING=UPPER XSF=1.0 YSP-0.5 ΕΙ; ΕΤΗ Τ᾽. 
MESH ZONE=PARENT NODES=0 ELEMENT=1 HIDDEN=DASHED BCODE=ALL 
* 


ADINA 
* 
END 


|) 


INPUT STRESS FN - THD3D4.STR 

INPUT ELEM. ANGLE FN = THD3D4.EL 

INPUT AREAS OF ELEMENTS FN =THD3D4.ARE 
OUTPUT FN =THD3D4.REL 

DUDEUI :THD3D€.REL 


- ہے ہے ہے‎ as aS ኤ سے سے سے سے سے سے ہے سے سے جے ہے‎ uo oo ወ x x መ መመመሙ ጠው ው ው = ው ው ው == መመ መመመፎመ Á un - سس كت كت‎ am am am am e =m =m am am am مس‎ =m am am ہت ہے ہے ہے کت‎ መው 


1 -0 1368704 00.660205 507477E*02 -0.669ET04  0.151E*01 :-09203E402 
EN-O00220E€04 ٧٥ 19292105 0,496ک‎ ١02 029ر‎ 04 ወ 2255٨7 د01‎ -0.381E*02 
SEO 2207/04 0. 199E+05  0,496ዮ+02 -0.2/9ET06 0.422 01 012 
ο ο 6:704د13‎ " 0.660EFE+05 ο. ΑΛ ΕΤΟΣ =0.669E+04 —=0.161E+01. 0.203E*02 
5... 220+00 7255105. 07 ]20E+t03F O 65BB5rtO4 :0.,267E*X02. 20  119ET02 
በከ 260+047 “0. 297E+05 0 ۰719+ 062-0245 ٢٢ ር 105817010: 105 0 د‎ 
EE ሇህ 1و‎ 0 ٢٢ l ۹05 ን 245 ۳۶039۰ >0 05401 02 
EE 2225300 R0. 725E 4050. 120ET03 —0.6BBETOS اح‎ 26 ٢102 -0.119E*02 
BEN θοε ο. 0.34 2E«-DBELSD.I55ETSO03490.5427E £06  .0.774E*02 -0.-334E*02 
NE 519-06. O,MJBEXONMEO.TISESOSCE0.422E405 J0:115E+02 0.387b*01 
ο... 115-06۷ ο ο τος 20 11 د‎ 0 20.627:2406. —0.115E*02 0. 387E+01 
ت06 ح60‎ 00363 5+ 051-0 ΤΕΒΕ 0 ሮስ 10. ስንን ያር እ ተሪ. 0.3 72 
۰۰۰۰۰۱766 0۷ 0۰70957005 =O. 985E+O02 0۰1822+08 0:4۔‎ 215+02 -0 4۹2 
νι Εν; ο ο 20.15 6 +03 0.792E*T04. 0.977E5E*01 -0.764E-0) 
٣ ους 0 200+ 057-0۰1564 6097۰ 07927۲068" ٣26079776۰017074 66+01 
Per Veco OSs ο 25 ος ο ΕΕΣ ۲0. 67 6۶047077 01ھ‎ ۶۰2 2 
pe ከበ”. 1 225 Oa ο. ες ο 0, 15285+02 >0.180ሮ+05 -0.172E+02 0.2 
M 2715104 ο ου 0S 0.1 522703 ٨-0-5337 4 0. 51BE*Ol1 0. 632bEf*02 
ከር οσοι» ο 0. 152051 ؛ 0 0.53 د‎ -0.5188 7 01 0.202 
2۰٠۰1 ενος 0 مه‎ 5 ٢ 169+037 >0 1B0EF04  0.172E+02 -0.362£4+072 
21 —0.342E+02 0.6 07205 0. 2575+02 -0.126E+04 -0.969E+01 1-2 
εντος δα 0.148BE+O2 0.94BE+04 0.214E+01 0.7122 
25. 0.234E«*04  0.45BE«*05 C.148E+02 ዐ. 958ጄ+05 >0. 214ጅ+01 -0.727E+01 
24 -0, 352=+02- 0, 6075+05 0. 2575+02 -0.126E+04 0. 969ጄ+03 -0,169ጅጄ+02 


25 -0,5405+0ሩ- 0. 2365+05 0©.169E+03 -0.125E+04 -0.836E+01 ۰377 +72 


1 


26 
24 
21 
2 
0د 
31 
9 


a æ a e ar or ہے سے تب جے جم ہے ہے سے کت سے‎ UND ہے جے‎ i x ‹» መ» መ- ہے‎ (UP x m m m q =ወ ብው ጠው ብው ዛው σα ج مس‎ መወ سے‎ መወ መ s s መው ሓመ 


16 
127) 
15 
5 
20 
21 


D 


D. 


.456E404 0 
.456E«04 0 
23895 4105 Ὁ 
.410E+04 0 
.447E404 0 
.447E+04 0 
.410E404 0 
51 

.660E+05 -0 
. 17160 O 
+٤٥ 
.6606+ O 
.725E+05 0 
27 +0 0 
27 2 ٰ 0 
۰72+ 
0 ο 
.254E404 0 
.A56rLE*04 0 
3 706, 
ASE OS RO 
2655405 FO 
-266E+05 0 
ISSE 05S O 
.444E+05 0 
ጋህ ጋ 05.170 
: 309E F05 

.444E+05 

.607E+05 <0. 


.950E«04 
. 950ጄ+0ዓ 
1.1 15 
, 5 
) 04 
.B585104 


26 5 


. 1 3 6+ 4 
..ء‎ 7+7 
, 9 
۸.6 09 
2 04 
+۶ 
. 511۳4 
. 32254 
.409E+04 
5435+05% 
۵ 
.409+04 
.168E«04 
. OS 
.116 964 
. 168+04 
+0 ۹ی‎ 


SUBBETOS 


188E+04 
13 2174 


34202 


.196 3 
>> 3 
.125E+04 
145015 
.1256E:E03 
3633ء‎ 


3ءء 


77002 و. 
.496E+02‏ 
ሳን UE‏ 
ο οπου.‏ 
-120E403‏ 
3 11. 
.119Ε03‏ 
.120Ε03‏ 
239 ۸. 
IIE TOS‏ 
۷ 
+7 ء, 
7ی 
1566+03- 
156E 0‏ 
2 
149E 03‏ . 
152153 
~LS2ZE+Ce‏ 
149E 03‏ . 


+32 


120 


ى٣٦‎ 
29003 
+66 OL 
9٭>7"ھھء‎ 
.600E+03 
.600E+03 


SIS SETS 


መር 45 
- 795E+04 
ك7‎ 
.669E+04 
.688E-*04 
.666E+04 
.666E+04 
.688E+04 
.542E+04 
.410E+04 
.410E+04 
:542ET04 
428ھ‎ 
.180E+04 
.180ጄ+ዕሩ 
.181ξ 104 
«180Ε504 
224 
223E C04 
.180E+04 


.126E+04 


1 ۷۹ 
+7 
و‎ 7 ES 
ο ΕΣ οὐ 
. 94 35+0ጊ 
"94 717 


۷ت 


. 203E+02 
LOU SEXUS 
379E F072 
20 SEXUS 
+2 
-LOSE+OZ 
ο ο 
-1195102 
.334E+02 
ب7٦1‎ 
GILE 1 
233453027 
“αυ τος 
ሳር ጋ OL 
6٦ 
.409E+02 
. 3627402 
«4305102 
.430E+02 
ΕΕ ος 


16 9ፔ+02 


.127 ٣ 
. 12753 08 
. 469E 02 
. 241E ٢" 
. SOE 0I 
0 7 


>2 ٦۷ 


16 -07 
.578E+01 
2ے‎ 
.161E+01 
.26 7/571۴ ۶۷۹ 
.lB4EdX M 
.184E«01 
.262Ë10_ 
.774E*02 
ء٣‎ 
.101E30 
7۶/۱۹٢٣٢ 
٭>ء‎ ھ٥۳‎ 
~114E+02 
. 111 10 
.421E*02 
.172bE4*02 
51/5091 
ء‎ 7۳۲ 
ء‎ ۳۶ 


θε... 


75۰ 
το 
0 
90: 
755 
or 


907 





22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
E 


ο -- - s ጩ‏ ہت ہت ہے ہے ہے جب ہو  —  — H s‏ ست ست پس مس مس s‏ ست مس مس ست مس مس پس مپس ست پس نس نس نس መመመ መመመ  ሙመሙመሙመሙመጩመጩመጠ‏ س ست کس کس سه s‏ سه መመ‏ کس ہی نس مس سب ټس پس — ہی — — ہی سے سی سے سی سی سی سے سے س س - سی سے s s‏ سے ہے ہد 


10 
11 
12 
JC: 
14 
15 
16 


Oe 


9 


D. 


σος 


0890 6 5 
70+1۲ ۶05 
NEO? Uo 
22955405 
.7 +94 
.870 +4 
πο τος 


. 265+05 


741E+04 
741E+04 


265505 


.000000 


. LOCE+01 
.993F+00 
.992E+00 


πο ο 


ምግ 
Ol 


0 
2100E-*01 
91 
ΟΡΟΙ 
B UDE-O01 
0ى‎ 7 
۹3002 3 
11 01 


.100E+01 


το. 


ο 
5075103 
39 2 
.44 0+4 
. 3.1۹ 
377+4 


. 238E+05 


410+4 


. 330+04 
. 330+04 


.410E+04 


πο 71 
1 G 
210052. Ol 
2100FE+01 
.906E«00 
3725100 
ο) ο. ος 
. 0 
.614E+00 
.133E+00 
413 3800 
26142400 
-996E+00 
"۹+ ۵ك‎ 
ο ας ο 


.996E+00 


0. 


۶0 


158+0 2 


τοσο» 
no ESOS 
۸۹03ء‎ 
.196E*03 
"00۶۹3 


(25+04 


ΕΠΟ 
¿J 26F105 
,. 3 
ISE 05 


7111 ፣)| 
۔‎ 1 
ο ο 
σος οι 
BEOOETOI 
1005701 
..71 
+٦ 
COE TOL 
210023400 
. 1001 
. 71 
ھ٦‎ 
"71 
.ء‎ 7۶01 


,7٦ 


191 


.265Ε04 
7 6 4 
6ء‎ 4 
.125E-*04 
.326E*04 
.326Ε04 
۹ءء‎ 7 
.193E-«04 
3+ 
"SABES OS 


.193E+04 


-100E- 01 
.۷[ 
111101 
5፡1) 11፡1] 
89۹3۹ 71 
+0 - 1[ 
[ΕΠΟ 
.100E+01 
"1005401 
DUOGUET01 
+٦٣ 
ο δν. 
... +17 
0 
+۷۷٦ 


ھ٦‎ 


1 
ο Ἔτοι 
ESET 
1502 
.276E+01 
۰27607 
. 468-72 
.21 +2 
2 7 
7.۔‎ 071 
2721 2 


20-7 
S095 4-01 
. 00701د‎ 
1 005+ 01 
PLOOE+01 
+1 00E +01 
2002401 
+7 
1007 
81۶ 
۸)7 
1 
۷ھ‎ 
7. 1.11] 
31ء‎ 


.100E+01 


‘187E O0 
06ء‎ 
. 71 
.86 0 
3 6۲٣ 
+7 
د‎ 2 
220295102 
1+02 
. 781E+01 
028902 


0 .۔ 
045r 000‏ 2 
0 4 
0 26+ 
sts oes OO‏ 
2769500 
76SE+00‏ . 
soe DoD‏ 
0.۔ 
.966E+00‏ 
፡ #66ር፣ፒ ተ00‏ 
0 ..-۔ 
ΘΕΣ ΕΟ‏ 
71+ 
1000 . 


7ء 


-663E+00 
+3 0 
. د‎ 0 
236535199 
»-661E+00 
.2B86E*00 
+26 +0 
+08 
«54 SEND 
1252-20 
ھ,‎ 0 
پ0‎ 
و‎ 0 
<999E+00 
“ 2 199 


22201 111) 


17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
26 
28 
29 
30 
31 
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210171 REL. 


— — — — جح سے نس پټ‎  — سے — سی‎ ጮመ።መሙር። መመ መመመ መመመ መ መ መ መ መመ መመ መመመ መመመ መመመመመመመመ መመ  መመ=  መ-ፎሬ- መመ ہے ہت بہت ہت ہت ہت‎ መመ መ=መመመጨመመመመ መመ سے‎ 


. 00 
ى7٤‎ 
10021 
.999E*00 
9975: 00 
.+0۵( 
. 999E+00 
. 9975100 
. 1005401 
. 1005101 
1 ንበ] 
ΠΕΠ 
100781 
. 1075 ۷ 
000 1 


"۷۱۳۱۰۷٦ 


0.960E+00 


ΠΠ 
.990E*00 
. 9905 +00 
.998E+00 
111111 
.100E*01 
. 10 08412 
ከከ ከ ከኮ 
.10 1 
1ος. οι 
LOGE Cl 
. 000E 00 
ΕΟ 
7ھ‎ 
sd ODE 111 


.100 +1 


0 0005100 


LOLAL STRUCTURE RELI P rr Y - 


=- መመ መ መሙ — < — — i ii ii መመ ሙመሙ፡ዱ ሙ መሙ = — مس‎ 
— — — — — — — — — ټس پا سا پا — — — — ټس — تی پس <= — — — — ټس‎ Z— ټی‎ = == = == = = <= = = 


07101 
ھ0‎ ٦ 
0.101 
0.10) ٦ 
0.10 1 
ο ۷٣ 
0.100E+01 
0.100E+01 
oc POOE FOT 
0 +7۲ 
0 ۶ء7٤‎ ٦ 
۷۹ھ "ء0‎ 
ο 1: 
ο 7 
O 1005 OL 


ο ο οἱ 


0 ۹۹ 


0.000E+00 
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ους Cm 
2110110120: 
111] 
.100E*01 
ους του 
.100EÉE+01 
.100E+01 
۔‎ 71 
SI QUESO 
. 72 
. 71 
LOO + OF 
.100FE:01 
"1095101 
. 1 00+ 1 


.100E*01 


0.000E+00 


+7 
. 10001 
۷ھ‎ 
[0 7 
800+7 
.+.. 7 
«1005101 
,. +7٦7 
.100E+01 
. 1 0892141 
7 LODE + On 
.. 07 
. 100E 071 
. 100 1 
LOSE EOL 


. 1001 


0.009E 00 


0ً 70ؤ‎ 
0) UD 
0.998EX90 
0.86 
0 ,18610)፻:1071 
0.997E*00 
0.7۲٦ 
0) 
OQ. OE 7 
0. SBE U0 
0.95300 
0.100E+01 
0.999B 079 
0. 985EX09 
Q 2985 ESI 


43 ۷۵ 


0) ۵ 


Q4 9:9 O0 
0. SESE 00 
0.988E*00 
Q.- S97 ESO 
0.99 7E ۷ 
0.996E+00 
0.9968400 
0.555. 
0.100E:01 
0.989 ۵ 
0.939. ۲ 
ο. οσα 0 
0.999E X00) 
0,985 ο 
0.3985 ٥ 


0. 999E 0 


0.000E+00 


እን ተክ 9 
TTTTTO 9 


TEE GOO 
Mad ce 





A. a 


N 


ك لال 43 --. 


0ل لال لا / لے 


Ye 7 

EE OT = ΡΟ 2 
az Aza -- 
Gz'mz- NIWAX 70 


ur MNSIUJJ 43 /ፎፎህገ6ህ-3 “34፡1 [1121፲ 3==> MEI ۰۴061 


ΕΟΤ Η 49): 


oc 


تئ۶ 


0٢11ا‏ - ۲11 683۸ 01ہ 





070 


25 


τττουυ 
434 =ሁ1 


143383 323 /43ህ/፡5 :3ሀ343 1 11 ος ο] 
ت٥6‎ Η αι ΟΕ 


3-3 “Ὁ 


ΠΗ! LAL) 

puc u- pP መን ን ን | 
رلاد 01 یں‎ r 
οφ οὐ - HIMHAX ገህክዘ፲ሟጧሟ፲፡0 


TINO 


٣٣٦٣‏ .| ےریہ 1160773 ٢ار‏ 3ت 


NOIGHZA WI-YNIGY 





UHIU: 
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κ ο ο ο A =o 1 ON ۶۳۱7 FETE 
x 


ለ ን 1141. .1ክብ EIGHT SYMMETRIC LAYER USING 3-D. 
አ 


FILEUNITS LIST-8 L06-7 ECHO-7 
CONTROL PLOTUNIT=CM 

COLORS BCODE=BLUE 

k 


TE BASE CREATE 
WORKSTATION DEVICE=0 
* 


meme ING * THM3D¢.IN; *AAEIGHT LAYERS SOLUTION 90,75 DEG ۰7 
ል 


MASTER IDOF=000111 
PRINTOUT VOLUME=MAXIMUM IPD=1 IPRIC=0 P=STRAINS 
THULE FORMATTED=YES FILE=60 


κ 

COORDINATES / ENTRIES NODE D Y 2 

1 -20.. -20. 0.3 

ο ٤٣٦ 

3 υπ 

4 -20. ο. ο 

5 0. DEO 

6 00 ο. 0.1 

7 -20. ZEE | 

8 0“ 75 

9 ο. ιο 
21 -20.  -20. ο. 
22 0. -20. 0. 
23 20.  -20. 0. 
24 -20. 0 
25 0. O On 
26 20. የ ο. 
27 -20. 2 MOS 
28 0. DEO 
29 ο. 20. 0. 
101 5 ο 
102 0. 17 
103 5. 17 
104 15. 2.68 0 
105 7.68 10 0 
106 15 10 0 
107 15 -10 0 


— — — ھی سے پس ہی ہے‎ — — — — — - -- — = — — — — πῶ — — — — — — — — — — — — — — — جم سپ‎ መመመ መመመ መመመ መ መ መ መ 


— — w — سو ټس‎ w — — መመ መመመ — سے‎ — — — — — — — — መ መ መጨመመል መ ጨሬጨመሬ መ።መጦጠ መ መ መመ መ መ-ሬ-መ መመ ټس‎ O 


MATERIAL 1 ORTHOTROPIC 181.E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 
۱ GAB-7.17E05 7.17E05 7.17E05 
MATERIAL 2 ORTHOTROPIC 53.8E05 17.9E05 17.9E05 0.083 0.083 0.25, 
GAB-8.6E05 8.6805 5 


۸ 
ο Μι CYLINDRIC X=0 Y=0 Z=0 PHI=0 THETA=-90. XSI-90. 
S I INATES / ENTRIES NODES R THETA XL 

10 105-9; 03 

11 OM D 

ገ 1011 ٣٦ 03 

15 3 nS on 01 

14 2 0 0:1 

5 200 ορ. UI 
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A 

LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 


LINE 
* 


۴۲۸۱۱ 1 THREEDSOCID 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 


CYLINDRIC 
δι ο, 
CYLINDI IIC 
CICINDRIC 
CYLINDRIG 
CYLINDRIC 
CYLINDRIC 
CYLINDRIC 
Cx DINDRIC 
CYLINDRIC 
۳سس ۶ت‎ 
CYLINDRIC 
CYLINDERS 
۲۳۲٣ 
CYLINDRIC 
"هي‎ 11. 


9 کی ف‎ QJ tJ eN 
-J CO tO OQ» QJ t ኑ፦” > 


16 


GSURFACE 7 8 28 27 
GSURFACE 8 9 29 28 
AXES-ORTHOTROPIC 1 
AXES -ORTHOTROPIC 
AXES-ORTHOTROPIC 
AXES-ORTHOTROPIC 
AXES-ORTHOTROPIC 
ASTRESSTABLE 1 1 
STRESSTABLE 1 21 
EDATA / ENTRIES EL 


2 
3 
4 
5 
z 


l 
ῬΈΕΙ ΤῸ 
29 
Z 
STEP 4 TO 
30 
S 
SIEP 4 I0 
311 
: 
SIEP 4 TO 


NODI=101 
NODI=101 
NODI=101 
NODI=101 
NODI=101 
ο ο 8:6٠٦ 


EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL-1 
EL-1 


ο. 
45. 
SO 


MIDNODES-1 
MIDNODES-1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES =] 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES =1 
MIDNODES=1 
MIDNODES=1 


270. 
315. 
350. 


1551 
1807 
225: 
270. 
317 
360. 


NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 


MATERIAL=1 TABLES 


EL1=1 
EL1=1 
EL1-1 
EL1=1 
EL1=1 
EL1-1 
፲,1,1 = 3 
EL1=1 


EL2=1 
EL2=1 
EZ 
EL2=1 
Ευ 2Ξ1 
EL2=1 
EL2=1 
EL2=1 


ርጋ ርጋ ርጋ ሮኃ ርኃጋር ο ο‏ تک ی ہے 


EL2=4 NO=9 
EL2-4 NO-9 
NODJ 2103 
NODJ=104 
NODJ=105 
NODJ=106 
NODJ=107 
27 


TABLE MAKES 


1 


1 
1 


ونم 


1 


1 
3 


GJ CJ 


Ü)‏ نب 
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NODK = 
NODK = 
NODK = 
NODK = 
NODK = 


HL HH 


EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 


1507 
107 
107 
102 
102 


NODES=27 
NODES =27 
NODES =27 
NODES =27 
NODES =27 
NODES =27 
NODES =27 
NODES =27 


ዶ 
LOADS ELEMENT 
9 
πο 
16 


7 


22 


2 


=3. 0504 


25. 0E04 


Lue 23 LINES / 8 28 
ΠΡ 23 SURFACES / 


E 28 'ٰ USES" 29258 


κκ 

κ 

BGROUP 2 THREEDSOLID MATERIAL=2 TABLES 

۳۰۰۰۰۱۱٣۶۱٦ 12-5 5 31 32 25 25 EL1=1 EL2=1l1 EL3=4 NODES=27 
"ከ 12 13 5 5 32 33 25 25 EL1=1 EL2=1 EL3=4 NODES=27 
۰۱ک‎ 13 14 5 5 33 34 25 25 ELI-1 EL2-1 EL3-4 NODES-27 
EXSBUME 14 15 5 5 34 35 25 25 EL1-1 EL2s1 EL3s4 NODES-27 
ΙΕ (۶ 5 5 35 36 25 25 EL1=1 EL2s1 EL3-z4 NODES-27 
۳۰۰۱٦٦ 16 175 5 36 37 25 25 ELI=1 EL2=] EL3=4 NODES=27 
۰.۰۰۱۶17 18 5 5 37 38 25 25 EL1s1 EL2=1 EL3=4 NODES=27 
EUM 101155303125 25 ELl1-1 EL2-1 EL3s4 NODES-27 
AXES-ORTHOTROPIC 1 NODI=101 NODJ=103 NODK=102 ` 
AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 
Peeo-ORTHOLTROPIC 5 NODI=101 NODJ=107 NODK=102 
ENLDEOSTABDE 112345678 21 26 27 


ENMDESSIABLE 1 21 


ολ / ENIRIES EL TABLE MAXES 


1 1 1 
9 9۶ Το 

27 1 1 

2 1 3 
SIEP 4 TO 

30 1 3 

3 1 3 
peer 4 TO 

3d 1 3 

4 1 1 
IEE 4 TO 

7 1 1 


x 

k 

EGZONE NAME-PARENT / 1 

EGZONE NAME=PATCH / 2 

| ۱٢۴5 HeEADING=-UPPER X5F=1.0 ¥SF=0.5 XFMAX=22.5 YFMAXK=17.0 

MESH ZONE=PARENT NODES=0 ELEMENT=1 HIDDEN=DASHED BCODE=ALL 

ο. HEADING=UPPER XSF-1.0 YSF-0.5 XFMAX-22.5 YFMAX-17.0 

MESH ZONE=PATCH NODES=0 ELEMENT=1 HIDDEN=DASHED BCODE=ALL 
* 

ADINA 

* 


END 


27 


پر پل ας‏ پل پلر پفر پل پل پفر پل لر پل پل پل پل پر لو لو لو لو پل لو لو لو لو لو پل ἡς‏ لو لو e‏ پل پل پل پل پل پا پل κ‏ پل پل پل پل پل پر پل لر پل k kk e‏ 


PROGRAM TRANS3 


k Ak k k k k k k k k k‏ لر لو پل Pe ας ας‏ لو لو لو لو پل k‏ پل پل پل پل پل پل پل re‏ پل k‏ پل پل P‏ پر لو hk p‏ لو 


8 
C 
C 
C 
C 
8 
C 
Ç 
C 
8 
C 
C 
C 
8 
C 
C 
C 
C 
C 
۸ 


τ. 


THIS PROGRAM TRANFORMS THE .,)1ክር.... አህ: re ebay: 
EXPRESSED IN THE GLOBAL STRUCTURE 00111 لال لال تت لل‎ 

THE PRINCIPAL MATERIAL COORDINATE. AND TREN, GhECULATES لال‎ 
RELIABILTY PER STRESS COMPONENT AND PER ELEMENT AND FINALLY 
EVALUATES THE TOIAL RELIABILITY. ALL TAT 9 ΙΙ ን 
SUBROUTINE 0/7 


NEL # OF WHOLE SRUCTURE ELEMENTS. 

NELP : f OF PARENT STRUCTURE ΕΠΗ 

SX,SY, , . SY2 = ህከ. ብ... IN ህህ ስህ የ AME > 767 

THETA > THE LAYUP ANGLESPER ELEVEN i. 

51,52,... 56 THE STRESS IN THE MATERIAL PRINCIPAL CO O D DE 
R1,R2,...,RO = THE RELIABILITY PERT SIKE S 7 
RT1,RT2,...RT6: THE TOTAL RELIABILTY OF STRUCTURE ٥٥٣٥ 
RE . THE ئ55‎ 168 0 ۷۲ 

P THE ۰۷ 

ALF THE SHAPE FUNCTION PARAMER FOR WEIBULL 


BEFORE RUNNING “THE PROGRAM VERIFY IN THE PROGRAM THE ANGLE THETA 
INPUT AND THE # OF ELEMENT ΝΕΟ 7 ο. THE MATERIAL PRU. 
GIVEN IN SUBROUTINE BETA. 


x hh hk hh hh A‏ پو پو پو لو پر پو پو ἡς ἡς‏ پو لو he e‏ لو پو پو پو ار ہار لو لو لو لو لو پو لو پو لر پو پو پر پل پر لو پر پو پو پو پو لو پل لو پر پو لو پو پر لر پلر پلر پل 


INTEGER NEL 

PARAMETER ΠΡ ος ۱۳532 ٢٣٢ ο ου ος 8. 

PARAMETER (ROWS=100,COLS=6, COLS1= 7) 

INTEGBRSCPDHEPAUONEIE) 

REAL*68 A( ROWS, COLS) ; BURROWS ΕΠΟ» 770 

REAL SX(NEL) ,SY(NEL) ,SZ(NEL) ,SXY(NEL) ,SXZ(NEL) , SYZ(NEL) 
REAL L ΠΕΙΛΙΠΕΤ, μποτ በ 1 ክህ. ከ αμ. 

CHARACTER^20 FNAME,FNOUT,FNAMEI,FNAME2 


WRITE (97%) 11٢ 1011 7117116 δα. 

DEAD (he, (are Ne 

WRITE (59%) ΠΕΙΣ 51185 ٢ 0070 

OPEN (UNITs15,FILESENAME,STATUSS OLD) 

WRITE {5:4}. Νε ΡΟ ΤΕ ΤΕ 
READ 79 ) 9 ۸ NAME? 

WRITE ) 5 ٭,‎ ( ۲۲8۷۲ ELEM ANGLE ΤΗΣ 
OPEN (UNIT=16 FILE=PNAMEZ ۷ OLD) 

NRIITE ο ο) ΕΕ ARES OE ٢٢٢٢٢٢٢٢ ٣ 
۲  ) (5) 7٣٤ 

OPEN (UNIT=17,FILE=FNAMEL,STATUS= OLE) 

ΠΤΙ OSES ENTER ΟΗΕ ٢٢٢٢-٢٢ 

RED CS CA ۹٣٢ 

OPEN (UNIT=13,FILESFNOUT 1210 1٣٢ ) 

WRITE (135%) “INPUT STRESS ٢ ጎጨ! 
WRITE (13,%) 'INPUT PLE ANGLE ENS ERNE 
WRITE ) ٭ رت3‎ (٣۰۴ ۲1160۰) ΠῚ ΠΡΙ 
WRITE (i374) ους εν ο ን እ ያከ 


^INPUT OP ANGLE (THETA) FORTEA ۶۳ 


DO 20 J=1,NEL 
READ )16:7( ٧,١ ۵:١ 3 4,۱۹,۱۱۶ ٧٢٢ 


& 112,113,114 ,115 115 ከ? 
EP 9) 7 )17ل‎ ۰0 ἘΠΕ, ٧٢ 
1 )11 ١۵) و(1‎ ١ 3 ΠΕΙ +) 0٣٥۷ 


WRITE SCS ۶۷۳۷۵ 
THETAP (T= ΠΕΕΑ PP Bone 
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20 CONTINUE 


ΙΤ; "OUTPUT :',FNOUT 
پل‎ 


n SL LORMATION OF STRESSES TO PRINCIPAL AXES STRESS 


πα نر وا دو ور‎ "542 SHY ,' SxEZ2 SY ANGLE 
οὐ ο I-l1.NEL 

ο ο (۱۱۲۴۳ 505 5S5X(-(1);5Y0(1), 526011), ,SZY(1);,52X2(1),SY2 (1) 
MEI) DPHEIAFCI)AI80.0 / PI 

|. ο G1) IEDLSAUIDTSYCI),S201) SZ Y CIJ,SAACIO,SYZ(I), 

& THETA(T) 


E- COS(THEITAF(IY) 
M=SIN(THETAF (I) ) 
ال‎ oe AM) IM SYULUIACHNSA2) 4 255RY(I)ALAM 
ህህ. ..፡፡1 κ σ ET οκ ΕΚ) — 2ASARYCL)ALAM 
= SZI) 
Seo lo ete oy lye la SAY IIA LAR? = Mano) 
ο ποτ ο ιο κ FE SALLIE 
ን፡፡ ۹ٔ ۹8 77701 
55 CONTINUE | 
ως ፡.፣ NEL’, Si Of foo 2 80 7 55 ,'94. 
POREO T=], NEL 
I I ο 1 ο لاکد لل‎ ο) Αι )دال در‎ 5٢ ۸۷٢1 1፡7 
80 CONTINUE 
* 


ZSUBHROUTINES FOR CALCULATING THE RELIABILITY 


ETT BETA(LA,B,NEL) 

Εμ σσ) ‘INPUT ALFA’ 
BEND (5,^) ALF 

υπ ο ο "ALDEA :',ALF 
ο ο ο δι DB NEL ያ 
CEL HEDI(A,B,NEL) 

 -. ከ1 1). 141) 


10 meme 70/1 له ۵7 0:6 ل۵‎ / 12 5600 =) 
ib BPCO IXOAIO, UCAT.ANO Ww o3X,800 —^) 
15 ጫው. ከ STIA TDHEIAU' I2, )s |) 

2 mem (ix, THETA Το لالب(‎ FO.) 

60 O i T (7/,4X,IG,0Ell.33) 

61 HU /,4A,IO,bEII.3,3X,F6.1) 


) 
) 


STOP 

END 
አተተተ+ተተ+4+4ተተ+ተ4+ተ+ተ+ተ+ተተ4+ተተ4+4+ተፅተተ4ተተተተተተተተተ+++ተተ++ተ+ተተ+ተተተተተተተተተ++ተተተተተ+ 
Jk A k k ہا اہر چا‎ oh oh Gh IA A پو‎ 

SUBROUTINE BEIACA,B,NELD) 
A C A چاو پل‎ GA AGI k A لو پل‎ Ph k K AI GR I I 

PARAMETER (ROWS=100,COLS=6) 

REAL^8 A(ROWS,COLS),B(ROWS,COLS) 

ΤΕ ο ο ο δις οσο ὅσα 
À 


En 7] 55 => STRENCIH,RAIIO CALCULATION 


١۸ IO) TEL SRI ۰6796۰69۳۰۰۰6۰ ۳ھ‎ 
EOF THE PARENT SIRUCTURE, 


٦ 


DOO 11. 
>= OIRENGTH PROEERTIES ت٣٦‎ 


€ 9 9 9 9 ده د‎ 9 9 9 G @ © GG 89 9  : 89 89 9 9 9 9 Q 9? ፃ ቀቅ ቅ ډو‎ 9? ፃ ፃ 9 9 9 9 ፃ ው% ፃ 9 9 4 8 


X(1) = 15, 054 

X(2) = 0.45ኞ4 

X(3) = X(2) 

X(6) = 0.6BE4 

X(4) = X(6) 

X(5) = X(6) 

X1C = -15.0E4 

0 ہ٦‎ 

5230 = 220 

IF (ACI) ۴ 
IP (AUR 2) 020008 E 
IF (AUI,3).]099 9000030 =e 


DO 80 J=1,COLS 
BC 1 pd) eS MO O 
80 CONTINUE 
WRITE (13,00) ΙΤ, ΒΤ) ΒΤ, ΖΒ, ΡΕ Τη; Ευ. 
90 CONTINUE 
xD, ۴08 IBE TALE 


N=NELP+1 
DO ٤ 
C STREIIGTIH ۲ ۶ ሆሆታ” ة٥‎ SE 


(6-6 6 ο سے نو مت‎ νο το ο τι. 
X(1) = 10. 344 
x02) ۷ٰ6 
7) 
X(6) = 0.414E4 
x(4) = X(6) 
ACS) = 17 
AICTE I FES 
አረር. 296990 
ASC = X2C 
Ε.Π. ۴ ب0‎ 
ΙΕ Γ.Σ ٢1 ee 
IF (ACLI3).bP 200 NOD 
DOSS 2-1, CULS 
B(I,J9 = ۶ٰ  + 15 


15 CONTINVE 


WRITE (13,60) I,B(1;1), Bil 2 3111 الک ٢ل لل‎ ٣٣ص‎ 
100 CONTINUE 


0 ΘΑ; ο 7ھ )076 1ظ‎ 9 b DUNS 0 
60 ۲۷۶۱۲۸ 4 8 ال‎ 

RETURN 

END 


Je لو پو لو پو لر لر پل‎ e ee پیر‎ e ہل‎ e e e he kh او‎ 
SUBROUTINE ALFA BNET 

Ke A A k k A k A K e eh پل پل پوپ پل ہل‎ K K K پل‎ k k KOK K OK K K K 
PARAMETER (ROWS=100,COLS=6) 
REAL^8 A(ROWS,COLS),B(ROWSOSSRSD 
REAL AREA(100),AMIN 

ü 


* PLY STRESS - STRENGTH,RATIO RI SED D لل‎ p 1780٦ 


WRITE (13,10) NEL, Μι yo veer Ve EN a: 
DO SO το VEL 
READG(Ə(17 58) 777 7 


190 


AMIN=400. j 
ITE CS 9S5) 1110, 0 
ΤΕ,» 55ο) ር!” ር  , 7۳ 
DO 80 J-1,COLS 
^ DNE UOS BODL) )WA*ALE 
AROE BCI ۰۰۰(۸ ٤٦ 
80 CONTINUE 
Ωω απ ο ο ΝΕ) Αι 2) δα ου) ο ο δα 5) δι) 
90 CONTINUE 
10 CEL IIR ALO, GOAT ARS) ΙΧ ΒΩ )) 
60 ο ከ {/,74X,16,6511.3) 
EIS ۲۰۰۱۷۷۷۰۰ |. ከከ... 12,204 2— "9219229 
RETURN 
END 


Fx o Ας ہار پار پار پار پار‎ k A K K A K k k k ہار پار ہا‎ A À A À A 
SUBHOUTINE RELICA,B,NELD) 

A hA A A Ao A AAA A 
PARAMETER (ROWS=100,COLS=6) 
REAL^8 A(ROWS,COLS),B(OROWS,COLS) 


τ STERS5S5 TENSOR RELIABILITY CALCULATION (WEIBULL) 


ከ! ከ)!” !1!1 ኖር) IRL TRL CRY TRG TRS URS 
DO 90 I-1,NEL 
POT BO J-i, COLS 
οι ο) τοι αι ο) 
80 CONTINUE 
ο ο 3 o0) ل1 0 2 د‎ 211 2) BCLOS)9$ BI 5) BCL, 5), Btl,4) 
90 CONTINUE 
ιο ο UIN ١11,6٨5, 45) 7 Ιλ ل١)‎ = )) 
60 ο ΡΕ ΔΝ ΤΘ SELI 3) 
RETURN 
END 
ἀλλά و پر پل‎ K Ó پل‎ k A A o Ae A e A Á K 


SUBROUTINE ELRE(B,NEL) 
لو‎ Ας Ας Ακ Ἂς Ας Ἂς Ας Ἂς Pe Rx x x XA پر لو لر بل لو لر لو بل‎ 
PARAMETER (ROWS=100,COLS=6,COLS1=7) 
REAL*8 B(ROWS,COLS) ,C(ROWS,COLS1),RT 
x 
AELEMENTS RELIABILITY CALCUL., AS WELL AS "FIBER" (R1) AND "MATRIX" 
ARELIABILITY CALCUL.& TOTAL STRUCTURE RELIABILITY CALCUL. 
ο مت‎ 0) ει, R1 9452/7 R37 “RE =’ RS”. RO’ 6’ RE’ 
Rr =1.0 
ο 9ی‎ 
DO 90 I=1,NEL 
C(I,COLS1) = 1.0 
DO 80 J-1,COLS 
CET IBI J) 
. ο ο οι CI COSI CT) 
LLI IO ۰ J) 
80 CONTINUE 
Pal - ከ11. 0151) 
TIE CC ICI CEILI CII CIL 4 CCI, SI, CCI OCC.) 
ΑΙ ο ο I C ١ 1 27 7 ο πο ο σας πο απο 
90 CONTINUE 
τήν. ''; 10; TOTAL REL ἘΠ, RIZ RTI RI6, ‘RIST, 
& 'RT4' ,'RT' 
ΠΕ ES CNEL ο Πριν ο) U COINEL 21 3) C(NEL+1 4), 


|! 


10 


65 
70 


C(NEL+1,5),C(NEL+1,6),RT 


WRIIE ΠΠ 
WRITE (13 7 OR I 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
RETURN 
END 


κ διοτι ۹ 4:7-٥٣ 

)/ o TEILT) 

(uU SEE) 

(4/713, TOTAL STRUCTURE RELI ADI ILI I ን 


185 


δικο τις ο. 


7 


ከ ከ! SIRESS FN = THM3D4.STR 

MEUT ELEM. ANGLE FN » THM3D4.EL 

ከ1 AREAS OF ELEMENTS FN =THM3D4.ARE 
OUTPUT FN =THM3D4.REL 

ን ከህ :TRM3D4. REL 


መመመ መ ሙ መመ መ መ መ መ መ መ መመ መ መመ መ መ መ መ መ መ መሙ سے ست ست مس سے سے‎ ሬሙ። መ መመ ሬመሙመ መመመ መመኤዴዱ- መ መቃ» መሙ ሙዴ- መጨ سے‎ ወው መዉ - መ መመ  ሙመ ጭ፡ጭ፡ መጨ سے‎ መመ سے سے مسل مس ست‎ መሙ መ ہے‎ m d πρ 


1 -0.277E+03 0.391E+05 0.187E+02 -0.235E+04 0. 909ጄ+00 0.622E+01 
pao. 279F+03 :0.209E*05 21600 0.236E*04 0,849ጅ+00 `0.604E+01 
3 0.279E+03 0. 209ጄ+05 0.ጊ360ጄ+02 0.236E+04 -0.849E+00 -0.604E+01 


۶۰٦3 0.391605.. 0. 1B7E+02 -0.235604 207909Et00 =0.622EŁ01‏ 0:2 وؤ 


ο ጋጌ BEt03 0.23291E+05 O0.248BE 101 =0.235E104 0:. 23299800 Ο.675Ε101 
νι ፡!|035 0.209EF05. 0.217870) 052356404 0.361E+00 0:681 

7 1.1.5 . 209E*05 ο. οι 24 -361E+00 -0.6B83E+01 

8 ος ο τα ο ος JC BREST 11 . دد24‎ 4 .399ጄ+00 -0.675ጄ+01 

9 ሪከ. ከኃ . 05 ን መን ሪጋ -552640 : 2 
10 ΟΠΕ د0‎ ھ٤۶‎ . 746E+00 +6 4 553+401 ۵ 
11 5 . 20385 . 746E+00 .235E+04 >9 ۱۸۰ ء-ء‎ 2 
12 .27BE*03 .391E+05 .653E-01 , 4 9 ۹ 2 
1: 52370603 SIPETO . 11-1 . 22322004 ΟΕ ΟΙ ህከ ትሪ 
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51 +1 >۰ 005+07 - ٣٦٣ SL OOE+ O1 ۵ھ آ ۱ھ‎ +1 0 
22 .., 90 00 በ 00.01 . 1005101 1005501 2 0 .994 +0 
53 ο ο οὐ 0ء‎ +7٦ ΠΕΙ 008 ٦ ۵ؤ 9و9‎ .994ጄ+00 
54 SUE 0I 1081ھ 1005101" 0787 100ء۰‎ .+۷۵ -9 +۷۵ 
55 100E OL .. 1 71+ 11114 01.1. ። 41401 Ετος .975F4100 
56 ΘΕΟ , 00237017 0 0 0+37 ΠΕΠ 3995ET00 ο 00 
51 ΟΠΕ 0 APODETOI 0 1005401 -۔‎ ۷٣ .. 71 ο Ετος ΗΠ 0٤۱ 
58 ΟΕ 1 .10 1237401 0.1 100E TOL ἔπ ΒΕΠ JIBE 00 
S9 ግ ከ1111 )۳ء‎ Or 91 , 1) +07: "Ε.Ο +00 . 2 75E+00 
60 “0٤0 ከስስ ፡ ፡ ህን ۰ء‎ OT 00.101 ο Εσοι 227 ን ٤۵ . 99۹٢0 
en E O0 -1006+01 T0. 100E+01 . 1005-01 . م100‎ 2 9 +٤ -994E+00 
62 211111 ΡΕ ο! 1010. 073 ,110 11 71ء‎ ×× SESTO +1 0 
63 ΠΡΈΠΟΙ ۰ 072787 .. 7٦1 ΠΕΙ +0 +0 
64 ۰۰۰۰ ህከስ ο 5777 HL OORT OL . 71 10ہ‎ .46 +0 
TOTAL REL. RIT RIZ RIS RIG RIS RT4 p 
ο το L OOO 0005100 0 00011700 0 0008700 0-000E*00  0-000E£00. 0:103E*00 
TOTAL ENNRUCTURE RELIABILITY = 0 1035: 00 


== Z= = Z= Z = Z =Z Z = — — — — — — — = s πι ہے ہے‎ ia ei _جسے تہب ےھ‎ — s O ےہ سے‎ 
πο w s s s s s — — s s ° s s s s s ہےھ سه سہھ نی سی سس‎ "em EP — = = بس سې سه سه سه‎ = ss 


Jo 


ττττττ J 





πι E 
τττουῦ 
ايل 31ؾ‎ 
F 
ae το 4 49 4 4 0 
۱ ው 33383 EE 
2 60 07 نو اګ‎ 
2 εως πρι. dUNIOSISO  UNIGU 


wee [D6 O5] 'ጋገኳዊዘህ ]ዘ434934 1343343410 “Sea. 1H9I3s« "NI ‘Oe 
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k ΕΤ Νο I N TENTE 36 ΠΠ; 
x 


ο. 4. IN, EIGHT SYMMETRIC LAYERS, REPAIR W/ DIFFER. LAYUP ANGLE. 
* 


FILEUNITS LIST=8 LOG=7 ECHO=7 
CONTROL PLOTUNIT=CM 

COLORS BCODE=BLUE 

* 

DATABASE CREATE 

WORKSTATION DEVICE=0 

* 


ING THA3D4.IN; ^^^EIGHT LAYERS, DIFFERENT PATCH ANGLE, [90,90] *X*^' 
* 


MASTER IDOF-000111 


PRINTOUT VOLUME=MAXIMUM IPD=1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED-YES FILE-60 
ደ 


COORDINATES / ENTRIES NODE X y 2 
1 کہ‎ 
2 በበ... በ 
3 200۶ ۶۰ 0-1 
4 -20. 0. 0.1 
5 0. REOR 
6 20. s 
1 -20. η RETE 
8 +0 7 
9 20. πη 

21 ንግ ይበ ED 

22 ጠጠ 9 2. 

۳ ھ7 23 

24 -20. ከ በይ 

29 (ju 017 0 

26 20k 177417 

2; -20. ο... 

28 0. ET 

29 20. 2ου ο 

101 5 G ο. 

102 0. ο... 

103 S. ο. 

104 15. 2.68 0 

105 7.68 10 0 

۱ 106 15 10 0 
107 15 210 0 

S 25 S 2 U ምም ሚመረ መጠይ 55552 5225252 2ک‎ 


ہہ ټی پټ پټ ټسټ — — — — — — — — — جب ہے ہے — — — — — — — — — — — — - - — — — — — — سے سے ہے سے  —‏ پټ ټسټ پټ ہے — — یس 


— መመመ መመ — سے — پس سے س‎  — — سے سے س س پس پس پټ ټسټ پس پس‎  — — — — جس‎ መመመ -ፉ መ -መ ቁወመፌመ — — — سے س‎ s s ټسټ‎ መ መሠ 


Eur ORTHOTROPIC 181.E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 
| کا‎ 0:2717 007۰۰117 ۴ ΕΤ. 


* 

SYSTEM 1 CYLINDRIC X=0 Y=0 Z=0 PHI=0 THETA=-90. XSI=90. 

COORDINATES / ENTRIES NODES R THETA XL 
10 0077 ος 
iJ 10. 45 0.1 
12 ο 5ο 0.1 
13 joe 125 0 
14 10. 180 ο] 
15 JO. 12 IE 0.1 
16 ο. απ 0.1 
1.7 10. 315 0.1 


Á 

LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
LINE 
x 


EGROUP 1 THREEDSOLID 


ΟΡ 07 
یں‎ +۳ 
CYLINDRIC 
CYLINDRIC 
CYLINDRIC 
CYLINDRIC 
ΜΠ 
CYLINDRE IC 
CILINDRIC 
CYLTNDRIC 
CYLINDRIC 
) ۶۶۳ 
CYLCINDRIC 
CYLINDRIC 
CYCINDRIC 
CXEDNDRIC 


EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
EL=1 
1 
EL=1 
EL-1 


MIDNODES-1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES-1 
MIDNODES-1 
MIDNODES=1 
MIDNODES=1 
MIDNODES=1 
MIDNODES-1 
MIDNODES=1 
MIDNODES=1 
MIDNODES-1 
MIDNODES-1 
MIDNODES-1 


Jen 


NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 
NC=A 


خش رک تچ رح خم oO‏ ته ርጋረ ን‏ 


P 


GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
۶ئ‎ ٦٣ 
۶ +06, 
አ! ۹ 1 
رد )۸۷۷ھ‎ 
ho OR OLR Or ۲ 
AXES-ORIBOTROP IG. 4 

ጋ 

2 


CD tO O^ با رم نا‎ ቴ +J 
ړل‎ OO (O (፲ነ ጋህ ኪን ኑ” کہ‎ 
ت۲‎ 
e 


BXESSODUSDESROPIC 
tO DRESS Bier ΤΙ 
5115513١ ٢ ٢ 
6 “+۹ ከዓ. 


1 
پت‎ ۵ 
29 
2 
STEP 4 τὸ 
30 
3 
STEP 4 TO 
31 
: 4 
31 To 
22 


"۶" TABEES 


EL1-1 
EL1=1 
EL1=1 
EL1=1 
EL1-21 
EL1=1 
EL1=1 
EL1=1 


EL2=1 
پاظ‎ 
EL2=1 
EL2=1 
EL221 
اط‎ 1 
EL2=1 
EL221 


EL2=4 NO=9 
EL2=4 NO=9 


NODI=101 NODJ=103 
NODI=101 NODJ=104 
NODI=101 NODJ=105 
NODI=101 NODJ=106 
NODI=101 NODJ=107 
3 (6S 16s 8 7 eae 


TABLE MAXES 


i} 


1 
1 


= 


1 


1 
3 


(ል) ሂኢ) 


UW‏ لم 
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NODK= 
NODK= 
NODK- 
NODK= 
NODK= 


2 


EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 
EL3=4 


1072 
1592 
TOZ 
2 
102 


NODES 227 
NODES=27 
NODES=27 
NODES 7 
NODES=27 
NODES=27 
NODES=27 
NODES =27 


LOADS ELEMENT 
9 
πο 
16 


* 


2 


22:9 01504 


46 


E523 LINES / 8 28 
ETE SURFACES / 7 828727 / 8 9 29 28 
* 


< 

۶۰۰۰٠٠٠٠٠٢۷٢ 2 THREEDSOLID MATERIAL=1 TABLES 

۰۰۰٦۴۶ 11 125 5 31 32 25 25 ELl1=1 EL2=1 EL3=4 NODES=27 
GVOLUME 12 13 5 5 32 33 25 25 ELl=l EL2=1 EL3=4 NODES=27 
[IE 13 14 5 5 33 34 25 25 EL1=1 EL2=1 EL3=4 NODES=27 
GVOLUME 14 15 5 5 34 35 25 25 لد للم‎ EL2-1 EL3-4 NODES-27 
۰۰۰۰۱٦۱515 l6 5 5 35 36 25 25 EL1=1@EL2=1 EL3=4 NODES=27 
۰۰٦۶۱٢ 16 17 5 5 36 37 25 25 ELl=1 EL2=1 EL3=4 NODES=27 
NE 17 18 5 5 37 38 25 25 EL1-1 EL2-1 EL3-4 NODES-27 
EDBUME 1J0 11 5 5 30 31 25 25 EL1=1 EL2=1 EL3=4 NODES=27 
۰۰٠٠۰0۶۲۱٢1۶0106 1 NODI=101 NODJ=103 NODK=102 
AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 
AXES-ORTHOTROPIC 5 NODI=101 NODJ=107 NODK=102 
 ΙΛΠΙΕ 11234567 Β 21 26 27 


ENHEDSTABLE 1 21 
ES ENIRIES EL TABLE ΜΆΧΕΣ 
J: l 1 
STEP 4 TO 
29 1 1 
2 1 l 
SIEP 4 TO 
30 l 1 
3 l ሺ 
DIEP 4 TO 
31 I 1 
4 l l 
EREP-4 TO 
Seb 1 1 


نه 
Ἂ‏ 


ο ο ΝΑΜΕΞΤΗΠΕΡΡΡΘΟΙΙΟὉ / 1 / 2 

ΓΙ... 1 IING=UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
DESHSZONESTHREEDSOLID NODES-s0 ELEMENT-1 HIDDEN-DASHED BCODE-ALL 
* 

ADINA 

ኣሓ 


END 
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INPUT STRESS FN = THA3D4.5TR 

INPUT ELEM. ANGLE FN = THA3D4.EL 

INPUT AREAS OF ELEMENTS FN =THA3D4.ARE 
OUTPUT FN =THA3D4.REL 

OUTPUT :THA3D4.REL 


NEL SX SI SZ SAY 54 SYA ANGLE 
1 -0.336E+02 0.371E+05 -0.470E+00 -0.197E+04 0.142E+01 0.316E+01 20.0 
2 0.690E+03 0.216E*05 -0.521E«00 0.295E+04 0.133E+01 0. 322E UM 15.0 
3 0.690E+03 Ο.216Ε105 -0.521E+00 0.295E+04 -0.133E+01 -0. 322E TOM Hor 
4 -0.336E402 .0.371E*05 -0.470E*00 -0.197E*04 -0.142E401 -0.31G6EX0m 90.0 
5 -0.722E+03 0.382E+05 0. 137E02 =U 2038104 -0, 12368701 CGC ZEAE 9۹9 
6 -0.759E+02 0.215E+05 Ο0.15ΒΕ102 0.279E+04 -0.666Et00 UO. SOE 5.0 
7 -0.759E+02 0.215E+05 Q0.,158E+02 0۰279۲+04" 0666+ 0٠ 5-0 
8 -0.722E+03 0.382E+05 0.137E+02 -0.203E+04 0 213871 2 p P DRE 2 0 
9 -0.530E+03 0.393E+05 ' 0.120EFE+402 20 2201704 0. 201800 .ሪ ያ ን Ρ EE 

10 0.46 88102 0.213E+05 0.137E*02 0.248E«*04 0.433E+01 0 1 P ٦ mss 

11 0.468E+02 0.213E+05 .0-137E*02 0. 2461104 -0. “551101 CO 1 oe 70 

12 -0.530E+03 0.393E05 0.1070 ο. 1285704 ο 8011-06 ٤٣٦ 90 

13 -0.1B9E+03 " 0.397E+*+051-0.200F102 0 241F ٣۴ 0۶۱٠۷۹۰۹ 50۴۰۲۷۷٣۷۰۶ 0۰۰ PD DE 900 

14  0.380E«03 0.213E+05 -0.166E+02 0.23) 0 Όρο. ος 750 

15 0.380E+03 0. 2325+05 -0.166E+02 0. 238804 UO πο οι 0 ٧٧۹٢۹ ٣ 5 8 

16 -0.189E«03 0.397E+05 -0.200E+02 -0.241E+04 -0.390E+01 0. 153ጅ+02 90 

17 0.166E+02 0.378E+05 -0.109E+02 -0.223E+04 -0.117E+01 -0.141E+01 90.0 

18 /0.627/E+03 .0.208Ef05 “ο Σε ο) 0 د0830‎ 0 11۴67١1 ο 3 ms Ὁ 

19 0.627E+03 0.208E+05 -0.882E+01 0 250681704  0.-172EX010 10-3715 0 T30 

20 0.166E+02 0.378E+05 -0.109E+02 -0.223E+04 0.117E+01 0.141E+01 2070 

21 -0.564E+03 0.377E+05 0.347E+02 -0.209E+04 0.981E+00 -0.526E+01 300 

22 0.411E+02 0.209E+05 0.344E+02 0.260E+04 0.167E+00 -0.618E+01 750 

23 0.4llE+02 0.209E+05 0.344E+02 0.260E+04 -0.167E+00 0.618E+01 50 

24 -0.564E+03 0.377E+05 0.347E+02 -0.209E+04 -0.981E+00 0. 526ጅ+03 90.0 

25 -0.293E+03 0.375E+05 -0.101E+03 -0.243E+04 -0.309E+01 -0.679E+01 90.0 
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26 
27 
28 
29 
30 
31 
32 


10 
11 
πο 
13 
14 
15 
16 
7 
18 
19 
20 
zl 
22 
23 


.146E*03 
"1 46F+03 
953+03 
8 9 9 ٤+ 0 
.664E+03 
.664E+03 
Ετος 
973 5+ دن‎ 
۰27+03 
9 275+03 
87000۴ 
8 74E+02 
ھ94“‎ 3ET02 
Ετος 
.874E+02 
E598ET02 
ΕΚ 
۷6 +02 
3698E402 
"i +2 
٣ 75+02 
"51 7E+02 
.465E+02 
.822bE*02 
ΕΕ ΕΤΟΣ 
ο 02 
۰22+02 
ου ΕΤΟΣ 
ΟΡ ሰር 


.130E-*03 


ΠΡ 
Ετος 
ο 5 
+۰ 5ء‎ 
2227E*05 
τες» 
116505 
95 
006 
+108 |) 
.308E+05 
ο 5 
221114 05 
.301F: 05 
፣ ጋ ትህ) 
። 4 1 ተሊ!) 
“+93 3 ٤ 
ο ος 
.300E405 
BOSSES 
፡.ሪሥ ۶ٰ 
ΕΕ 
“505 
ء٣۳‎ 
8+36 05 
πο δν» 
30 2EËE+65 
SUSE US 
+5 


.309E+05 


ο 
=. 
ο 


Ore 


το 


84+72 
844+2 
MOTE+Cs 


140E+03 


270 ل1 
1122E POB‏ 
.140E+03‏ 
2085102 
...1 5211 
ο ۳۶‏ 
20852 . 
.748E+02‏ 
7 61 : 
-613E+02‏ 
.748E+02‏ 
ን ከ 9 ELO‏ 
ማር. 72‏ 
Bg ERO‏ 
ο ۷‏ 
-410E+02‏ 
ΠΕΠ.‏ 
ጋሪ‏ 


2119 τα 


552 


ን 2 
00 2 
. ص۵2‎ 
+2 
"+0 ۶ ۹۳2 


2 


ος 


0. 


το 


πο. 


۲ 


0s 


= OF, 


το. 


187E+04 


187E+04 


.24 3+04 


ሪጋ 11.5) 


ራን ۶۹ 
. 265E*04 
.252E*04 
9 3 
ΠΕΠ 


9. 3 


538E+03 
690E+03 
690E+03 


690E+03 


SOSOEST03 
.775E+03 
ን 114 
Ag SEES 
+ ۷۲۹ 
SUES 
.-۔‎ 03 
+133 
۸ - ۹ 
+7 ት 
,, +59 
+6 ۶ 
ግ ከ ٰ 3 
.668E+03 
-GOBETO3 


+96 3 


1 5 


=O. 
zc 
ο 
ο. 
29 
ED 
0. 
Oe 
Dr 
0. 
20- 
ος 
0. 
ο. 
ED 


zi 


. 71 
.457E*01 
5055101 
.267E*01 
POZE On 
-502E+01 
. 2675101 
ከሪ! 13.1 
+9 76 
. 0 
¿2235 1 


EBD ELO 
3255100 
329F+00 
5025401 
3060 
124E+01 
124E+01 
230 OT 
150E+01 
2306 11) 
Ba EPOD 
150E401 
178E+01 
B09E+00 
ο ٤ 


GE + Ou 


.724E-*00 
.211E-«00 
. 2111 0 


. 4287701 
+8٦ 
+ 1 
.487E+01 
"2955101 
ΘΕ 
-487E+01 
.212E+01 
. ት] 
-۔‎ +1 
'-20  - ٤ 
ο آ7‎ 
396 ۷ 
+ 71 
. 3025101 
+917 
767ھ‎ 
AS ETOT 
Οδ το 
ΕΠΙ 
5202. 101 
ΕΤ 
٦ 
ΙΕ 102 
+۸۷٦ 
6ھ‎ ۷ 
+0 2 
"l و9‎ 
2 


.116E+02 


Loe 
5 
90" 
20; 
τς. 
Vis 
90. 
80. 
21). 
96. 
30. 
950 
90: 
go 
90. 
90. 
90. 
90. 
90. 
90. 
5ο, 
90. 
90; 
SIUE 
90 
30. 
50. 
"ο. 
50. 
20: 


24 
25 
26 
27 
28 
29 
30 
31 
32 


— جح سے عو‎ s ہھ ہے ہو‎ መመመ መመ መመ መጨ پس نټ‎ መ መ መ መ መ መመ መ سے سے پټ ټسټ نس‎ መሙ መ መ ہت‎ መመ መሙ = መ= መ= ጨመ መመ መ መ መ መ cmn ےہ سعت‎ = መ መ መ ጠው መሙ መመ መ መ ጩመ ወሙ ጠው መ سے‎ 


10 
11 


JUS 


14 
15 
16 
17 
18 
19 


4ሪር 2 
.162F103 
.6 03 
επ. 
. 1620 3 
. 1631103 
+, 0S 
.166+03 
1635305 


ο τρις» 
<. 2166105 
217 0S 
.37 5 
5ءء‎ 
ΠΕΡ 
21۹5 
. 382205 
765. 
:. 6 ο 
OE 
392 5 
397E 05 
211FE+05 
.211 5 
7۶ 
ΡΕΞ 
.207E+05 


. 207 E405 


=O. 


ο Εις” 
. 31505 
31585 
15۳د‎ 
1۹11245 
.314E+05 
ο. Εν” 
. 14+5 


.314E+05 


.336E+02 
.612E+03 
δε ο. 


τ» 


72 03 


۶ 
~249E +02‏ 
2 258د ر/. 
۲ 3102707 
2220702 
επ 5‏ 
Εν.‏ 
5 ھ. 
596+03 . 
9(9 . 
75 ھ+, 
۶ ۶+ 
73ء 


.730Ε--03 


ος 
(on 


O. 


0 


ο 
- 0. 
ο 
UE 


UE 


374+302 
563E+02 
456E+02 


~456E+02 


5626702 
605E+02 
699E+02 
699E+02 


δο 2 


-470E+00 
۵ 
0 
.470E+00 
1 1 2 
. 7۶ 
ھ7‎ ۷ 
۸+72 
1121020102. 
17E ον 
1 27ÉE+02 
-12067072 
- 200 E102 
.166E4*02 
. 166+32 
- 200E+02 
- LOSE+02 
882301 


:.7 
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. 6688702 
«Θ20Ε 905 
-620E+03 
.620E+03 
-620E 03 
ΕΠΕ 
.589E5+03 
. 589E«03 
. 5B9E«*03 


~197E+04 
. 266E+04 
. 266+4 
۶ھ(‎ 
.2035+04 
. 298+4 
. 298+04 
2004 
.226+04 
.316E+04 
. 316E«04 
. 228+ 4 
.241E+04 
.316E+04 
316+4 
.241E+04 
2 ۶ 
መጋ ተህ 


. 289+ 4 


0۵ وی 
77141 . 
7+17 . 
٢1‏ 11 . 
٦۷ھ‏ 
.297E+01‏ 
leer OF‏ 
1141511« 


ΗΠ 


:.۸۹ 7 
. 7ج د‎ 701 
د.‎ 45511 
23155 7٦ 
. 284E+01 
.470E+01 
.420Ε01 
„. 284E Q1 
.4160ጄ+02 
. 1485102 
. 1 4 و35‎ 2 
1162 
72ء‎ 
1505202 
.150E«02 
. 1586402 
.141E+01 
-403E+01 


-403E+01 


0.125E+02 
0-۔‎ 617+71 
-0.48B5E+01 
0.4B5E+01 
0.617E+01 
0.6137 
0.337E+01 
πο 8٦ 
-0.613E3101 


-0.142E+01 
-0.451E+00 
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